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SPECTROSCOPIC STUDY SMALL SAMPLES 
HOLLOW CATHODE DISCHARGE! 


Abstract 


Operating characteristics the hollow cathode discharge have been investi- 
gated qualitatively. sample less than ugm. element may used 
study the hyperfine structure the principal lines exposure few 
hours. The useful lifetime the sample the discharge about hr., and 
80% the sample can recovered after use. 


Introduction 


The cooled hollow cathode type discharge has proved very useful light 
source for the study the hyperfine structure line spectra (4). general, 
few milligrams sample material are required. However, Schiiler and Schmidt 
(3) found that, with microgram thallium liquid air cooled hollow 
cathode, the hyperfine structure the green line (5350A) could photo- 
graphed. McNally, Harrison, and Rowe (1) were able detect gm. 
fluorine analysis problem. view the interest the nuclear spins and 
radioactive isotopes, was felt that further investigation the 
practicability using small amounts sample the hollow cathode was 
desirable. 

Apparatus 

The experiments described made use hollow cathode discharge 
tube conjunction with gas circulating and purification system, 
power supply energize the discharge, and (iv) quartz spectrograph for the 
examination the light emitted the discharge. 


The spectrograph used was Hilger equipped with quartz optics. 
quartz lens with focal length cm. was used focus the discharge the 
slit the spectrograph. Eastman IIIF emulsions were used record 


the spectra. 
The power supply and gas circulating system were conventional design. 


Hollow Cathode Design 
The discharge tube (cf. Fig. consisted water cooled anode and cathode 
separated Pyrex glass spacer sealed place with red wax. Brass 


Contribution from the Department Physics, University Alberta, Edmonton, 
This research was grant from the National Research Council Canada, and 
one the authors (F.H.F.) was holder Bursary the National Research Council 


Canada (1949-50). 
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was the principal material used the construction, while the actual electrodes 
were aluminum inserts, which could accurately positioned the anode and 
cathode sections. order simplify the diagram, soldered joints have not 
been shown. 


Noteworthy features the discharge tube are: 
(1). The use aluminum electrodes. desirable that the electrodes 
coated with some refractory material decrease the spectrum the electrode 
material. the case iron this accomplished artifically oxidizing the 
surfaces contact with the discharge. Aluminum, course, requires special 
treatment, since the surfaces are oxidized immediately contact with air. 


Fic. cathode mount. anode; cathode; gas inlet; water jackets; gas 
outlet; tapered plug for withdrawal cathode; quartz window; pyrex spacer. 


(2). The use tapered plug coated with vacuum grease provides for the rapid 
and convenient replacing loaded cathodes the cathode section. This 
design was possible since water was used for cooling rather than liquid air. 
The use the latter produces difficulties when joints are sealed using high 
vacuum greases and waxes, since these crystallize low temperature leading 
leaks the apparatus. 

(3). Placing the inlet tube for the circulating gas near the quartz window 
produces stream gas through the anode toward the cathode, which tends 
prevent sputtered cathode material from reaching the window and ren- 
dering opaque. 

(4). The aluminum cathode finally adopted consisted aluminum tube 
mm. inside diameter, threaded receive aluminum insert, and provided 
with six radial holes 0.8 mm. diameter set around the circumference about 
mm. from the front face. These holes were connected groove which 
joined opening leading the rear the cathode. The circulating gas 
entered the front face the cathode and left means these holes. This 
design produces (2) gas near the front face the cathode, which 
tends prevent particles from the cathode load from reaching the anode, 
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being deposited there, and being lost.to the discharge. The cathode insert, 
which formed the back surface the hollow cathode, was the shape 
truncated cone, the flat tip being mm. diameter. was this tip that the 
load material was placed. Earlier models the cathode were constructed 
one piece, the discharge region being simply hole drilled aluminum rod 
the desired depth. The characteristics the two types are discussed below. 


Experimental 
Loading Samples 


When small amounts sample material are used hollow cathode, 
becomes important examine the effects different locations the sample. 
Using the one piece cathode the type used Schiiler and Gollnow (2), 
was found that the sample was most effectively excited when placed the 
wall near the front. Loading was carried out using small hypodermic needle 
deposit known small amounts solution. Difficulty was experienced 
placing the solution the desired location and then maintaining there 
during evaporation. was convenient,, therefore, use the cathode insert 
described above. When the sample was placed the tip and evaporated 
slowly the use the insert provided similar discharge characteristics those 
cathodes loaded near the front. The gas ‘‘window”’ formed the circulating 
gas front the cathode opening did not appear appreciably affected 
when the distance between the tip the insert and the front the cathode 
was about 1.0 cm. When this distance was mm., the lines the load material 
were intense initially, but disappeared almost entirely after about min. 
operation, indicating that the effectiveness the gas had been 
greatly reduced. distance 1.5 cm. was found produce discharge almost 
completely lacking lines the load material. Subsequent discussion will 
therefore deal only with tip-to-front separation 1.0 cm. 


Operating Conditions 

order avoid the anomalous fall cathode potential, estimated 
Tolansky (4) take place when the current density exceeds ma. per 
the discharge current used this research was ma. Smaller currents were 
found produce discharges whose light emission was lacking the load lines. 


The intensity the light emitted the hollow cathode was found 
increase markedly the spacing was reduced. spacing 
less than mm. was therefore used for the bulk the work. The pressure 
the circulating gas was adjusted give maximum light intensity all trials. 
Variation the rate gas circulation within the limits available with the 
diffusion pump did not appreciably affect the discharge. 


Discharge Characteristics 


The assessment discharge characteristics this work was limited the 
intensity the load lines with respect lines the circulating gas, and 
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the length time over which the load lines could observed. Materials used 
load the cathode were sodium chloride, magnesium chloride, copper sulphate, 
zinc chloride, and silver nitrate. each case the load applied was gm. 
the metal. was found that the more intense lines all these elements had 
intensity comparable the helium line 2829A after five hours operation, 
providing that the discharge had operated continuously. Magnesium and 
copper lines were still detectable one hour exposure after nine hours 
previous operation. observing the load lines after resumption the dis- 
charge, was found that their intensity was greatly reduced. Apparently 
large fraction the sample condenses the less active portions the system 
when the discharge stopped and not released when restarted. 


Since equal masses the metals were used all cases, the numbers atoms 
available for excitation were greater for the lighter elements than for the 
heavier ones. For instance there were about 2.5 times many magnesium 
atoms available copper atoms. Yet the copper spectrum was general the 
more intense. Examination the insert tips after use indicated that the copper 
remaining had been reduced metallic form, appearing thin film the 
tip. The sodium and magnesium tips were somewhat whitish indicating the 
presence salt some kind. From the absence any observable film like 
that salt oxide, the silver and zinc tips probably also had metallic films 
after use. 


appears, therefore, that the intensity effects are related the electro- 
positiveness the load metal with respect aluminum. The latter electro- 
positive with respect copper, silver, and zinc. Hence its presence tends 
reduce salts these metallic form. the other hand, aluminum electro- 
negative with respect sodium and magnesium, and these elements remain 
salts perhaps oxides. Thus enhance the intensity the sample lines, 
conditions should arranged that the sample element not chemically 
combined. 


Recovery Sample 


Using load radioactive material, was found that after three hours’ 
operation with discharge current ma. about 30% the sample was lost 
from the insert. Anode and cathode each showed about 10% the initial 
count, while nothing above background could located the rest the 
circulating system. Thus, should possible recover 80% the 
load material from the anode, cathode, and insert, even after hr. 
operation. 


Suitability for Interferometry 


The principal lines the elements with which the cathode was loaded gave 
satisfactory blackenings the photographic plates hr. exposure when 
interferometer was used. Experience with the Fabry—Perot interferometer 
(mirrors aluminized 80% reflectivity) indicates that, when interferometric 
studies are made, this exposure time should increased between 
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and hr. This time the order the useful life the sample the light 
source, and could reduced the use more sensitive photographic plates 
than those used these experiments. 


References 
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THE EFFECT ION BOMBARDMENT THE EMISSION 
FROM OXIDE COATED 


REDHEAD 


Abstract 


Some preliminary measurements the decay emission caused ion 
bombardment oxide coated cathodes are presented. shown that cathode 
life increases with emission density ion bombardment prevented. 


Introduction 


The constant evaporation emitting material from the surface oxide 
cathode results ultimately failure the cathode emission, thus ending the 
useful life the tube. Frequently cathodes lose their emission long before all 
the active material the cathode coating becomes exhausted; this second 
phenomenon termed cathode and caused foreign atoms 
ions entering the cathode coating and altering its emitting properties. 
experimental study the changes emission caused bombardment the 
cathode surface positive ions formed normal tube operation has been 
made, and the subject this paper. 


may created within vacuum tubes three main processes: 


(a) collision the electron stream with residual gases the 
tube. 


(b) Electron bombardment contaminated electrodes (‘‘anode 
the energy the impinging electrons excess the heat for- 
mation the contaminating material, this material will decompose 
form ions (3). 

(c) When the cathode operated under temperature limited conditions, 
oxygen evolved neutral atoms molecules (1). This oxygen may 
ionized electron collision and the positive oxygen ions will return 
the cathode. 


Processes and are the most important. 


known that oxygen ions are the chief cause ion poisoning, pressures 
oxygen being sufficient decrease the emission several powers 
ten (6). Oxygen ions entering the cathode coating are presumed combine 
with the free barium the coating and thus reduce the density the barium 
impurity centers. The saturated emission density from the cathode given 
the familiar equation 


/9 = Ae 


Manuscript received February 28, 1951. 
Contribution from the Radio and Electrical Engineering Division, National Research 
Council, Issued N.R.C. No. 2484. 
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where constant, the reflection coefficient, the temperature, and 
the density impurity centers. The surface work function constant and 
the energy interval only slightly dependent m». Equation (1) shows 
how the emission affected ion bombardment the assumption that 
causes reduction the impurity center concentration, Other factors 
besides the change impurity center density are operative causing this form 
emission decay. 


Oxide cathodes are most sensitive poisoning when they are operated under 
temperature limited conditions. This probably due the evolution oxygen 
from the surface the cathode, mentioned (c) above. fortunate, since 
tubes are normally operated space-charge limited, that under space-charge 
limited conditions poisoning occurs slowly. However under pulsed conditions 
tubes are frequently operated temperature limited and the effect ion bom- 
bardment clearly apparent when using long pulses. 


When emission can recovered from poisoned cathode operating for 
short period increased cathode temperature, removing the source 
poisoning material, the poisoning called reversible. this case the cathode 
recovers ejecting the ions and the emission finally re- 
turns the value had before poisoning started. 


Emission failure tubes electronic equipment due either loss 
active material, mentioned the first paragraph, form irreversible 
poisoning which closely associated with reversible poisoning. Reversible 
poisoning occurs continuously during the life oxide cathodes, the cathode 
being poisoned and reactivated simultaneously. postulated that there 
progressive failure minute areas the cathode surface reactivate, thus 
the effective emitting area the cathode decreases with time. This causes 
the commonly observed phenomenon decreasing anode current and mutual 
conductance. 


operating tubes potentials less than the lowest ionization potential 
any residual gases the envelope, and also below the breakdown potentials 
any material deposited the electrodes, possible ensure that ions 
are formed. The effect positive-ion bombardment emission has been 
studied operating two tubes the same envelope (to ensure equal residual 
gas pressures), one tube being operated with anode voltage below the lowest 
ionization potential, and the other higher potential. Thus the first tube 
shows the decay emission due causes other than ion bombardment and 
the second diode shows the decay due all causes. 


Long Time Effects 


-2.1. Temperature Limited Conditions 


The decay emission under these conditions quite rapid, decrease 
50% the initial emission may occur within few hours (5). 
The decay temperature limited emission was first observed using com- 
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mercial tubes. Double triodes with separate cathodes were employed ensure 
that the pressure residual gases would identical the two tubes tested 
under different conditions. Batches tubes were taken, one half the 
double triode was diode connected connecting the grid and anode, while the 
other half was triode connected. The diode section was operated the 
anode, the triode section was operated 300 The triode grid was held 
space potential. The temperature the cathodes was such that both sections 
were temperature limited, and thus the initial anode current through both 
diode and triode was substantially the same. The only difference operating 
conditions the two cases was that the cathode the triode section was 
subject ion bombardment, while the diode cathode was not. 


ANODE CURRENT 


TIME (HOURS) 


Fic. Variation temperature limited anode current (in arbitrary units) with 
tially Tubes were diode connected and Tubes were triode connected. 


The results operating batch tubes under these conditions are shown 
Fig. which representative all the runs under these conditions. From 
time zero hr. Tubes were operated diodes and Tubes triodes. 
Over this portion the curve the diode emission remained sensibly constant 
while the triode emission (i.e., the emission from the cathodes subject 
positive-ion bombardment) decreased rapidly first and finally approached 
steady state. time hr. the connections the two batches tubes were 
reversed and Tubes were triode connected, while Tubes were diode con- 
nected. The decay Tubes similar the earlier decay Tubes Tubes 
show slow reactivation which continued for about hr. before reaching 
steady state; this steady state value was only 75% the initial value emis- 
sion. The initial emission could regained operating the tubes in- 
creased cathode temperature for few minutes without drawing current, 
showing that the cathode poisoning for this relatively short operating time 
was entirely reversible. 
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The conditions for the test shown Fig. were 
290 v., 

Cathode temperature 


Numerous types cathode core materials and coating material have been 
tested this manner. They were all found give very similar results. 


2.2. Space-charge Limited Conditions 

Oxide cathodes operated space-charge limited are considerably less sensitive 
poisoning ion bombardment than when operated temperature limited. 
This may due part the fact that oxygen evolved space-charge 
limited cathode. 


NORMALIZED SATURATION CURRENT 


LEGEND 
DIODE CONNECTED 
TRIODE CONNECTED 


TIME (HOURS 


Fic. saturated emission under space-charge limited operating conditions. 


batch double triodes was connected the previous experiments. 
The diode connected section was operated 1.2 the anode and the triode 
200 v.; the grid voltage the triodes was adjusted make the triode anode 
current initially equal the diode current. The cathodes were operated 
approximately 1050°K. Saturated emission was measured intervals lower- 
ing the cathode temperature about 750°K. and measuring the anode current. 
The emission from the diode section remained fairly constant for least 6000 
hr., while the emission from the triode section fell until reached steady 
value about 4000 hr. Fig. typical results obtained several tests. 
The normalized saturation current refers the current zero time unity. 


The presence high percentage silicon the cathode core material 
causes rapid decay emission owing the formation high resistance 


All cathode temperatures quoted were measured with electronic temperature bridge operating 
the retarding field region, having accuracy better than (2, 4). 
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interface between the cathode core and the coating material. Fig. illustrates 
such case. this example the core contained 0.25% silicon and the tubes 
were operated low current density (10 ma. per both these factors 
being conducive the growth high resistance interface. Both triode and 
diode connected sections are equally affected the growth the interface. 


LEGEND 
DIODE CONNECTED 
TRIODE CONNECTED 


0.7 


0.6 


0.5 


NORMALIZED SATURATION 


TIME (HOURS 10°) 


Fic. Decay saturated emission under space-charge limited operating conditions when 
interface formed. 


STANDARD ERROR OF THE MEAN 


LEGEND 
DIODE CONNECTED 
TRIODE CONNECTED 


Time (HouRS 10°) 


Fic. with time the standard error the mean the anode current batch 
tubes. 


Fig. shows the variation space-charge limited anode current amongst 
group tubes expressed the standard error the mean. Initially there 
was greater dispersion amongst the diode connected sections than amongst 
the triode sections. Within the first 500 hr. the dispersion amongst the diode 
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connected batch decreased rapidly while the triode dispersion increased. After 
1000 hr. both dispersions increased the same rate. 


Effect Emission Density 


usually assumed that increasing the emission density from oxide 
cathode will result reduced cathode life. However, experiments have shown 
that cathode life increases with increasing emission density provided that the 
anode voltage below ionization potential (7). 


Special experimental tubes were used these tests containing two diodes 
one envelope. The cathode sleeve Grade nickel (0.05% silicon) was 


LEGEND 
—+—+— HIGH EMISSION DENSITY 
-O-— LOW EMISSION DENSITY 


NORMALIZED SATURATION CURRENT 


ec 
n 


Fic. Variation with time the saturated emission from two diodes operated different 
space-charge limited current densities. 


coated thickness 1/20 mm. Two anodes different diameter (1.9 mm. 
and 4.5 mm. respectively) were used that different space-charge limited 
currents could drawn the two anodes which were maintained the same 
potential. The temperatures the cathodes the two diodes were measured 
under operating conditions with optical pyrometer and were found equal 
within the accuracy measurement. 


Fig. shows the normalized saturation current function time for one 
batch tubes. measure the saturated emission was obtained depressing 
the cathode temperature. The data Fig. were obtained from batch 
tubes with anode voltage and cathode temperature 950°K. The 
initial space-charge limited current densities drawn the large and small 
anodes were ma. per cm.? and 625 ma. per cm.? respectively. will seen 
that the saturated emission the diode operated low emission density de- 
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creased rapidly steady value about 35% the initial emission. The 
saturated emission from the diode with high emission density varied con- 
siderably the first 5000 hr. and finally remained constant about 110% 
its initial value. These general results were confirmed number batches 
tubes operated the same manner. 


has been found that the difference rates change emission greater 
low cathode temperatures. the cathode temperature was increased the 
difference the two rates decay decreased until 1000°K. the difference 
small. these high cathode temperatures the rate decay emission the 
low emission density cathode was always little greater than the rate decay 
for the high emission density cathode. 


postulated that increasing the emission density causes increase the 
electrolytic production free barium within the coating and thus enhances 
the emission. With voltages above ionization potential the increased current 
density creates large ion current which poisons the cathode and obscures any 
effect due increased electrolytic action. However low anode voltages 
ion poisoning occurs, and cathode life found increase with emission density. 


low cathode temperatures the production free barium thermal pro- 
cesses reduced and thus the electrolytic production free barium becomes 
the predominant process. This evidenced the large difference rates 
emission decay between low and high emission density cathodes low cathode 
temperatures. the cathode temperature increased, and the electrolytic 
process becomes less important relation the thermal process, the difference 
between the rates decay diminishes. 


Conclusions 


has been shown that positive-ion bombardment one the chief causes 

the decay emission from oxide coated cathodes. satisfactory fit the 
observed decay curves can obtained the basis the Wilson semiconductor 
theory assuming that the poisoning process due only reduction the 
impurity center concentration. Thus appears that this form emission 
decay due two more simultaneous processes. suggested that the 
decay emission caused reduction the number impurity centers 
the cathode coating and sputtering the cathode surface the relatively 
heavy ions. suggested that the final steady value emission represents 
equilibrium between the poisoning and reactivation processes. 


has also been shown that the life oxide cathode not diminished 
when the emission density increased factor much provided 
that ionization occurs within the tube. 


Investigations are continuing clarify the mechanism ion bombardment 
poisoning and will published later. was felt that brief description the 
preliminary experimental results would useful this time. 
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THE MACROMOSAIC STRUCTURE TIN 
SINGLE 


TEGHTSOONIAN? AND BRUCE 


Abstract 


Single crystals high purity tin grown from the melt modified Bridgman 
method are shown partitioned into bands, striations. difference 
orientation exists between the striations such that the lattice one striation 
may brought into coincidence with the lattice its neighbor pure 
rotation about axis parallel the specimen axis. The properties the 
striations have been shown dependent both the rate growth and the 
crystallographic orientation relative the direction heat flow. tentative 
explanation the origin this effect advanced terms the formation 
edge type dislocations from the condensation vacant lattice sites. 


Introduction 


Mosaic structures occurring metal other crystals may divided 
broadly into two main categories: micro- macromosaics. Micromosaics are 
those which occur order magnitude about one micron, while macro- 
mosaics are usually visually observable order magnitude. 


Considerable attention has been paid the past micromosaic structure. 
This has its origin the theoretical work Darwin (7) who suggested the 
existence within crystal mosaic block structure slightly differing 
orientation from block block order that observed and calculated X-ray 
reflection intensities should closer agreement. Since Darwin’s work much 
experimental and speculative work has been done which summarized 
paper Buerger (1), who has also pointed out the existence what 
called This was term applied crystal that par- 
titioned into branches each lineage ‘being descendent from 
common nucleus, and having slightly different orientation from its neigh- 
boring lineage. 


This paper describes some work that has been done type lattice im- 
perfection that produces mosaic structure the macroscopic type. Because 
their macroscopic character, these mosaics are readily observable and their 
behavior and properties may examined relatively simple techniques. 


careful study these mosaics has been undertaken single crystals 
high purity tin. They not, however, occur solely tin. They have been 
observed lead, zinc, copper, and silver single crystals well. Because its 
low melting point and relative ease handling, tin was chosen for more de- 
tailed study. 

received February 1951. 
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1935, Greninger (9) reported the results work that had done the 
lineage structure copper single crystals. grew single crystals copper 
the Bridgman method and examined the fine structure obtained the 
X-ray back-reflection Laue spots. From this fine structure was able show 
that the single crystal was partitioned into lineages between which there 
existed difference orientation from five minutes two degrees. also 
showed that these lineages had predominantly common [110] direction and 
that they were separated predominantly common plane. 


The present work tin somewhat resembles that Greninger, but differs 
from three rather important respects: First, the method growing the 
single crystal specimens such that crystals any desired orientation may 
grown. Even though Greninger’s method yielded high percentage good 
single crystals, did not permit control orientation. Second, the method 
growth yields single crystal specimens which the lineage structure 
visually observable. This provides very powerful additional way examining 
the behavior the structure under various conditions growth. Third, the 
rate growth may controlled and determined leisure after the crystal 


completely grown. 
Method and Materials 


Thus far, only high purity tin has been used for the experiments. This tin 
labelled supplied Messrs. Johnson, Matthey London, Eng- 
land. Manufacturer’s spectroscopic analysis shows this 99.987% pure 
with the principal detectable impurities being Pb, Bi, Fe, Sb, and Cu. 


The method used for growing single crystals tin was developed Chal- 
mers (4). The tin melted graphite container, called boat, which 
generally rectangular cross section. Typical dimensions are cm. long, 
cm. wide, and 1/2 cm. deep. This boat clamped horizontal position, 
extending into small electric furnace made few turns Chromel wire. 
The furnace mounted carriage which may driven any desired 
speed along tracks parallel the axis the boat means geared drive. 


seed crystal known orientation placed the boat the end remote 
from the furnace and junction made between this seed and the molten tin. 
The crystal then made grow moving the furnace along the tracks 
direction away from the seed. This will cause the interface separating solid 
from liquid advance into the liquid and eventually all the liquid tin will have 
solidified. all the operations are done carefully, the liquid will solidify into 
single crystal the same orientation, relative the axis the specimen, 
that the seed. 


Obtaining seed crystal start with presents great difficulties. 
this, some tin melted into boat previously described. Then, using glass 
rod, the liquid tin drawn out from one end fine and narrow point 
possible. slowly moving the furnace direction away from the pointed 
end, the molten tin made solidify. Conditions such this are very favor- 
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able for the nucleation, the fine point, single crystal which will grow 
the exclusion other crystals. 


Having obtained single seed crystal this manner, crystals any desired 
orientation may grown orienting the seed appropriate direction 
when making the junction with the melt. Single crystals tin have been suc- 
cessfully grown this method rates high cm. per min. rates 
greater than cm. per min. there strong tendency for stray crystals 
form. 


valuable feature this technique that periodically the liquid end the 
system may agitated. This produces ripple which freezes the instan- 
taneous position the solid—liquid interface. From these ripples the shape 
the interface and the rate growth the crystal may determined 
leisure. 


tin, the structure, which has been called becomes apparent 
after short etch solution 10% ferric chloride, 10% hydrochloric acid, 
and 80% water. 


Experiments 


has been pointed out that the technique used should allow growth 
crystals any desired orientation. practice, however, the success with 
which single crystals may grown is, considerable extent, dependent upon 
the orientation. has been found that the direction growth that 
which the direction longitudinal heat flow parallel the [110] direc- 
tion. Accordingly, most the preliminary experiments were done crystals 
which the [001] direction was perpendicular the top surface the crystal 
and the [110] direction was parallel the specimen axis. 


When crystals this orientation are etched, they show the top free sur- 
face series fine bands running down the length the specimen parallel 
the specimen axis; see Fig. These bands are quite strikingly visible when 
the surface examined with light appropriate angle, owing the differ- 
ence reflectivity each band. The widths the bands, striations, vary 
considerably, but are generally the range from 1/4 mm. mm. 


all but few cases, striations are not visible the bottom surface the 
crystal. possible explanation for this given below. 


Fic. Schematic diagram crystal boat. 


Usually crystals are grown from seeds that are narrower than the crystals 
themselves, smooth junction being made means graphite insert 
shown Fig. 
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PLATE 


Top surface single crystal showing striations. axis perpendicular plane 


paper, and axes 45° specimen axis. 

Fic. Fig. with 15° rotation crossing the specimen are one 
minute intervals. 

Fic. Fig. with 26° rotation about axis. 
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Fics. 10. X-ray pictures along one pair new striations taken mm. intervals. 
Note gradual approach the spots each multiple. 
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Under these conditions, has been found that the striations present the 
seed continue through into the crystal. This not possible for that portion 
the crystal that line with the graphite insert. this region, new striations 
form independently from the seed. can seen Fig. these new striations 
generally start some distance from the sloping edge the graphite insert. 


Crystal Structure Striations 


ordinary back-reflection X-ray Laue picture from single crystal consists 
single spots arranged hyperbolas. All these spots lying one hyperbola 
are due reflections from crystallographic planes one zone. 


the back-reflection pictures taken from the surface crystal showing 
striations, the pattern shows each reflection multiple instead single. 
The multiplicity each reflection depends the width the striations and 
the size the X-ray beam the surface the crystal. Thus, the multiplicity 
found equal the number striations covered the X-ray beam. 
Fig. shows the nature this multiplicity. 


All the spots given multiple are due reflections from planes the same 
(hkl) but each tilted slightly different angle the direction the incident 
X-ray beam. From the regular manner which the multiplicity develops for 
all reflections, apparent that there orientation difference between the 
lattices neighboring striations consisting pure rotation about common 
direction. This rotation, which would bring the lattice one striation into 
coincidence with that its neighbor, may measured using Greninger’s (10) 
chart, and has been found the range from approximately 1/2° 
4°. 


Crystals which are grown with the [110] direction parallel the specimen 
axis and the [001] direction perpendicular the top surface always show stri- 
ations parallel the specimen axis. all these specimens, back-reflection 
ray pictures show that the difference orientation between neighboring stri- 
ations rotation about the [110] direction axis. 


Series X-ray pictures have been taken across the width crystal over 
tions. These have shown that the orientation any point this region always 
lies somewhere range bounded 6°. This indicates that the rotation 
the lattice going from one striation the next not always the same 
direction. Rather, there are about many rotations one direction there 
are the other. 


This has been shown more precisely taking series X-ray pictures 
intervals across the width the specimen equal the diameter the X-ray 
spot. Then, superimposing neighboring pictures, set continuous meas- 
urements rotations the lattice from striation striation may made, 
both magnitude and direction. These results show that the sum all the 
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rotations one direction very nearly equal the sum all the rotations 
the other direction. 


Striation Distribution 


Two experiments have given conclusive evidence that the striations extend 
throughout the whole volume the crystal. 


First, crystal approximately mm. thickness was started growing the 
usual manner. With the interface advancing slowly into the liquid, 
the boat was suddenly tipped that the portion still liquid ran away. This 
decanting process leaves exposed what was essentially the inter- 
face the time decanting. Several back-reflection X-ray pictures taken 
random positions this decanted interface all showed multiple spot patterns, 
even though visible striation structure was present. The direction the 
multiple spots again showed difference orientation consisting rotation 
about the [110] direction. 


Second, crystal which showed very marked striation pattern its upper 
surface and virtually striations its lower surface was subjected very 
strong electrolytic etch. The electrolyte used was 10% solution hydro- 
chloric acid water and the current density was about 0.1 amp. per sq. cm. The 
over-all effect was reduce the crystal about one-half its original thick- 
ness. all stages etching, visual examination showed striation pattern 
each newly exposed surface. Back-reflection X-ray pictures from the final 
surface showed multiple spot patterns similar those from the original un- 
etched surface. 


Thus, there can doubt that the imperfections that give rise stria- 
tions extend throughout the whole volume the crystal. 


INTERFACE 


LiQuID 


Fic. 11. Schematic diagram shape interface. 


interesting note that, after etching, the newly exposed surfaces the 
bottom side the crystal also showed striation where none had been seen 
before. This would suggest that the interface quite sharply 
curved just the bottom the boat, where contact made with the graphite 
owing the higher thermal conductivity the graphite than the tin. 
Previous shown that, when the rate growth becomes suffi- 
ciently high, the trace the interface the surface becomes curved con- 
cavely towards the liquid. When this was the case, the striations longer 
formed direction parallel the specimen axis but instead they converged 
showing tendency grow perpendicularly the interface. Thus, the inter- 
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face just the bottom the boat was curved, the striations would tend 
develop perpendicularly this interface and would grow into the crystal 
rather than parallel the surface. (See Fig. 11.) The etching, having removed 
the thin layer where this had happened, revealed the striations layers further 
up, which had developed the normal way. 


Striation Width 


For crystals the previously described orientation, had been observed 
qualitatively that the width the striations increased the rate growth 
decreased. attempt was made correlate striation width and growth rate 
quantitative way. number factors tend make this correlation 
rather difficult: 

The difficulty growing crystals under exactly controlled 
conditions, 

The fact that striations any one crystal vary considerably width, 

The difficulty simultaneously observing and measuring the striations. 


The first these may reduced maintaining thermal conditions 
steady possible and waiting for equilibrium established between the 
solid and liquid before starting the timed growth. The rate growth then 
very nearly the same the rate motion the furnace. overcome the 
second difficulty, the average large number measurements taken. 
Also, instead listing average widths for single growth rates, the average 
width for range growth rates given. measure the striations, was 
eventually found best photograph the crystal and make measurements 
enlarged print. 


TABLE 

Rate interval, Average width, 
mm. per min. mm. 
0-1.0 1.20 
0.5-1.5 1.05 
1.0-2.0 0.85 
1.5-2.5 0.75 
2.0-3.0 0.60 
2.5-3.5 0.50 
3.0-4.0 0.45 
3.5-4.5 0.40 
0.35 
9.0-10.0 0.25 


Measurements were carried out crystals grown rates from mm. 
per min. down less than mm. per min. Table shows the results. This 
table shows trend towards constant average width the rate growth 
increases. This not strictly true. Owing the curvature the 
interface the higher rates growth mentioned before, the striations rapidly 


converge and become fine that they cease possess individual identity and 


any width measurement becomes impossible. 
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Stability 


The stability the imperfections giving rise the striations has been 
tested subjecting specimen successive 120 hr. and 240 hr. anneals 
227°C. (melting point tin 232°C.). From visual examination and from 
ray pictures taken from the same positions the surface before and after the 
anneal, has been shown that there change either boundary positions 
orientation difference. This shows that the structure the imperfection 
extremely stable. 


Orientation Effect 


All the experiments described far have been done with crystals oriented 
with the [110] direction parallel the specimen. axis, and the [001] direction 
perpendicular the top surface. order determine the effect the crystal- 
lographic orientation the formation striations, was necessary grow 
crystals with the direction heat flow parallel other crystallographic 


directions. 
109 


10 O10 


00 


Fic. Schematic diagram showing orientation crystal simplified stereographic pro- 
jection. Arrow indicates direction striations. 


With the [001] direction still perpendicular the top surface the specimen, 
can specify angle between the [110] direction and the specimen axis 
(direction longitudinal heat flow) (See Fig. 12). Thus, for the crystals pre- 


Crystals were grown which the angle was progressively increased 
steps degrees from 45°. For each orientation, crystals 
were grown various rates. all cases, the general results are the same. 
low rates growth, the striations formed are virtually parallel the specimen 
axis. the rate growth increased, the striations form progressively 
greater angles the specimen axis. Figs. and show this effect. The maximum 
angle which the striations make with the specimen axis increases the angle 
increases but not rapidly. Thus, there tendency high rates growth 
for the striations parallel the [110] direction. these high growth 
rates, the angle between striations and the [110] direction increases the 
angle between the [110] direction and the specimen axis increases. The re- 
sults are shown graphically Fig. 13. 


Here, the results are plotted the angle between the striations and the 
specimen axis. number features should noted these curves. 
the angle increased, the rate growth which necessary produce 
constant angle formation the striations increased. 
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the angle increased, the maximum angle between striations and the 
specimen axis also increased. 


Back-reflection X-ray pictures indicate that the difference orientation 
between neighboring striations always rotation about axis parallel 
the specimen axis, regardless the angle between the striations and the 
specimen axis. 


z 

x< 

3.0 
4.8 g=45° 
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° 4 8 12 16 20 24 28 32 36 40 
RATE per min. 


Fic. 13. relating direction striations and rate growth for different orientations. 


Formation Striations 


X-ray survey new striations that form the crystal independently 
those the seed has given the following results. series back-reflection 
pictures was taken along one pair striations line with the sloping shoulder. 
These striations were not visible for distance about one centimeter from 
the sloping shoulder. The X-ray pictures showed that, the region where the 
striations were well developed, there was fixed difference orientation 
between them. the shoulder was approached and the striations became 
less distinct and disappeared, the orientation difference decreased until the 
two spots could scarcely resolved. This shown Figs. 10. 


Although this effect has not yet been fully investigated, leads some 
promising possibilities. Qualitatively, has been observed that the lower the 
rate growth, the greater the distance from the sloping shoulder the 
point where the new striations are visually observable. For sufficiently high 
rates the new striations form right the sloping edge. Thus, would seem 
that, when new striations must form, the case when the geometry the 
arrangement makes the crystal wider than the seed from which grown, 
they require incubation period develop their maximum extent 
dictated the conditions growth. 
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Discussion 


Although the experimental data available are perhaps insufficient warrant 
very exact theoretical analysis, they suggest possible explanation the 
origin this phenomenon. Throughout this paper, continual reference has 
been made the production single crystals whereas the resulting product 
has always consisted members differing orientation. However, inasmuch 
the orientations all the members specimen are such that they have 
common direction and may brought into coincidence with that any other 
member rotation very few degrees, seems justifiable refer the 
specimen single crystal containing irregularities that result what could 
described limited fiber structure. 


The geometrically most convenient way that has been suggested describing 
the boundary that separates regions differing orientation means 
array dislocations. Burgers (2) has discussed the dislocation arrays that will 
describe the boundary that separates regions between which any arbitrary 
difference orientation exists. the present case striations, the boundary 
could described array edge type dislocations suggested 
Burgers alternative array edge dislocations suggested Shockley and 
Read (14). should borne mind that the dislocation models that have 
been proposed and analyzed are strictly applicable only simple cubic lattices, 
whereas tin based rather complex body centered tetragonal lattice. 
spite this, felt that, first approximation least, consideration 
striation boundaries terms simple edge type dislocation arrays would 
very profitable. 


this picture striation boundaries has any basis fact, there still 
remains the question the reason for the formation dislocations the first 
place. 


When solid crystal heated, will acquire concentration vacant lattice 
sites, which, according Frank (8) will the order its vapor concen- 
tration for simple molecular crystals and considerably larger for metallic 
crystals. This concentration will greatest the melting point. Calculations 
Huntington and Seitz (11) have shown that for copper the melting point 
this concentration may high 10°. the metal crystal cooled 
from its melting point, the concentration vacancies will decrease. Nabarro 
(12) has suggested mechanism for the elimination vacant lattice sites 
which the number vacancies reduced their migration and escape 
from the surface the crystal. However, the mobility vacant lattice sites 
considered too small for such mechanism operative any but 
the very smallest perfect crystals. 


second mechanism that has been suggested Frank (8) and Seitz (13) 
seems much more probable. the temperature lowered from the 
melting point, the vacancies, lacking other suitable sites eliminate them- 
selves, will form into aggregates. Both Frank (8) and Crussard (6) have indi- 


| 


TEGHTSOONIAN AND CHALMERS: MACROMOSAIC STRUCTURE TIN 379 


cated that the most stable form for such aggregate the shape flat 
disk. Condensation vacant lattice sites into flat disks would continue 
that, the temperature decreases, the disks would grow size. When 
critical size reached, collapse the lattice across the disks occurs, producing 
general ring dislocation. 


This would the situation the whole the solid metal specimen were 
kept the same temperature during cooling. the method growing crystals 
from the melt used this work, the circumstances are quite different. Here, 
the solid contact with the melt the interface the melting point and 
temperature gradient exists from this interface back into the solid direction 
parallel the specimen axis. The magnitude this gradient along with the 
rate advance the interface will govern the manner which the vacant 
lattice sites will condense and the subsequent collapse the disks. The effect 
this temperature gradient that the disks vacant lattice sites would not 
all grow and collapse throughout the whole volume the specimen the same 
time, which would the situation gradient existed. Instead, the solid 
immediately behind the interface, the disks would formed and exist such. 
small distance beyond the point where, virtue the temperature gra- 
dient, the temperature considerably lower, the disks would have collapsed 
forming pairs edge type dislocations. the interface advances, more vacant 
lattice sites would condense the disks the interface, and those parts 
the disks which have now fallen lower temperature will collapse, con- 
tinuing the propagation the pairs edge type dislocations. 


the relatively high temperatures immediately behind the interface where 
these dislocation pairs are formed, the dislocations would fairly mobile. 
Cottrell has pointed out (5), these dislocations would have tendency form 
into transition surfaces consisting arrays dislocations one sign. Cahn 
(3) has utilized this mechanism describing the single 
crystal specimens that have undergone small amounts plastic deformation 
followed annealing. This would result just what observed; i.e., instead 
single crystal, one obtains specimen which divided into parallel regions 
with small differences orientation between neighboring regions. Also, because 
the equality production dislocations opposite sign, one would expect, 
fact the case, that taken over sufficiently large number bands, the 
rotations one direction would balance those the opposite direction. 


this general basis, the effects changing the rate growth the crystal 
and the orientation the crystal relative the direction longitudinal heat 
flow may considered. 


The size which the disks vacant lattice sites will grow depends both 
the temperature and the time for which given temperature maintained. 
The temperature will control the degree mobility the vacancies that per- 
mits their agglomeration into disks, and the time limits the interval which 
the mobility the vacancies sufficiently high make aggregation possible. 
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Therefore, with given temperature gradient the solid, the smaller the 
period time allowed for the aggregation vacancies, the smaller will the 
size the vacancy disks. Consequently, higher rate advance the inter- 
face will result collapse the lattice occurring across vacancy disks that are 
smaller than for lower rates growth. 


addition, the higher rates growth correspond shorter times which 
the dislocations can assemble into transition surfaces. Thus, these transition 
surfaces will necessarily closer together. This produces narrower striations 
which agreement with observation. 


When crystal grown from seed which this mechanism has already 
produced striated structure, that part the crystal that line with the 
seed will grow with propagation the striated structure into the crystal. 
that part the crystal not line with the seed, new structure must 
initiated. the basis the mechanism described, one would expect find 
incubation period for the establishment stable structure with maximum 
orientation difference between striations. This too agreement with ob- 
servation. 


If, very probable, there exist preferred crystallographic planes into which 
vacancies condense, would expected that the actual manner their con- 
densation would depend the orientation these planes relative the 
direction longitudinal heat flow. Thus, the formation striations would 
direction that compromise between crystallographic and thermal con- 
ditions. has been shown that, when the [110] direction parallel the 
longitudinal heat flow, the striations form parallel this direction. However, 
when the direction makes some angle with the specimen axis, the direction 
the striations compromise between the direction longitudinal heat 
flow and the [110] direction, being parallel the former for low rates growth 
and tending closer the latter high rates. This would indicate that the 
vacant lattice sites prefer condense planes closely associated with the 
{110} planes. 


conclusion, should again emphasized that the theory suggested 
the preceding section is, best, only very qualitative. does, however, seem 
borne out fairly well the experimental evidence that thus far 
available, and itself completely self-consistent. 


Summary 


single crystals tin grown from the melt, type lattice imperfection 
causes the crystal partitioned into striations. These striations, which are 
present throughout the whole volume the crystal, are such that the lattice 
one striation disoriented from that its neighbor rotation from 1/4° 
about the specimen axis. statistical fluctuation exists the measured 
properties these striations, such their widths, their direction for- 
mation. spite this, regular variations these properties have been noted 
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functions the rate growth and the orientation the crystal relative 
the direction longitudinal heat flow. 


mechanism proposed account for the formation the striations. 
This based the condensation vacant lattice sites into disks, and the 
subsequent collapse these disks, forming pairs edge dislocations. These 
dislocations than migrate into transition surfaces stable arrays dis- 
locations one sign, producing the striation boundaries. 


Observations indicate that the {110} planes are preferred for the formation 
vacancy disks. 
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DENDRITIC GROWTH 


WEINBERG? AND BRUCE 


Abstract 


Dendritic growth large single crystals lead revealed decanting the 
remaining liquid any required stage the process solidification. 
direction dendritic growth was determined function crystallographic 
orientation and direction solidification. was found that all the dendritic 
rows single crystal were not quite parallel, owing crystal substructure. 
explanation advanced for the presence the substructure terms den- 
dritic growth based their observed properties. measurements 
during solidification showed that the region the liquid where dendrites were 
formed was supercooled several degrees. supercooling was observed 
the portions not showing dendritic growth. 


Introduction 


Although recent years considerable attention has been paid the nature 
crystal growth, very little has been reported concerning metal crystals grown 
from the melt, particularly for the high rates growth which dendrites have 
been observed. This probably due the opacity the melt which prevents 
the direct observation crystal growth, possible with some organic materials, 
and the higher temperatures involved with metals. The authors have there- 
fore examined the dendritic growth large single crystals lead, and will 
present some the observations. attempt made this paper present 
any theory explain the formation dendrites during growth, nor will exist- 
ing theories discussed. 


Cross sections dendritic structures have been observed for many years 
polished and etched sections castings and are familiar most metal- 
lurgists. The structure consists skeletons solid metal, formed during 
solidification, which grow into the liquid either ahead the general 
liquid interface from solid particle the liquid. For cubic metals, their 
shape usually described that Fig. adapted from the diagram 
Carpenter and Robertson (1). This structure has recently been questioned 
Northcott (3) his examination copper alloy castings. 


From the figure, the dendrite seen consist primary horizontal stalk 
from which small secondary stalks grow periodically the horizontal and 
vertical planes. From these secondary stalks, smaller tertiary branches grow 
parallel and perpendicular the primary stalk and on. believed that, 
under proper growth conditions, fine primary stalk initially forms from the 
solid. Secondary branches then grow the stalk becomes longer and thicker, 
followed tertiary branches, the entire skeleton becoming progressively 
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larger and thicker growth proceeds. The direction growth the primary 
stalk appears depend the crystallographic orientation the parent metal, 
Northcott and Thomas (4) having shown for copper alloys that dendritic mark- 
ings polished surface were outcrops {100} planes. 


Since expected that dendritic growth may explained terms 
atomic arrangements the solid adjacent the liquid, terms heat 
flow, both which are markedly influenced impurities, only very high 
purity (99.999%) lead was used. was also apparent that, there was regu- 
larity the dendritic pattern, should more readily observable with large 
crystals, and therefore relatively large crystals were grown and examined. 


Fic. Dendrite skeleton. 


Papapetrou (5) and others have advanced explanation the dendritic 
growth nonmetallic crystals from solutions terms diffusion phenomena, 
the essential condition being the occurrence supersaturation the liquid. 
analogous manner, they believe that dendritic growth for metallic crys- 
tals may associated with supercooling the melt. 


series cooling curves was therefore obtained for various stages growth 
the crystals order obtain further experimental evidence relating super- 
cooling dendritic growth. 


Experimental Procedure 
Method Growth 


Because the high purity the metal used, the crystals could not pol- 
ished and etched reveal their dendritic structure, since this method implies 
the presence second constituent impurity the metal along the dendrite 
boundaries. Further, confines all observations the plane polish, showing 
only section three dimensional dendrite. alternative method was 
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therefore adopted. This consists pouring off the remaining liquid the de- 
sired stage growth. The solid that has formed remains the boat and can 
examined leisure. 


The crystals were grown the manner described Chalmers (2) which 
the ‘‘seed single crystal the desired orientation, placed 
horizontal graphite boat containing some polycrystalline metal. This metal, 
and part the seed adjacent it, then melted small furnace sur- 
rounding the boat, and the liquid agitated make good contact between 
liquid and seed. The furnace then slowly moved away from the seed along 
tracks parallel the boat. The resulting crystal generally single crystal 
having the same orientation the seed and grown rate governed the 
rate which the furnace moved during solidification the metal. seed 
crystal different orientation grown above, but the original seed 
tilted and rotated the required amount with respect the boat axis. The 
boat consists strip thick graphite which trough 3.4 cm. wide 
and 1.0 cm. deep has been milled, leaving side walls 0.2 cm. thick. The furnace 
mounted wheels and consists dozen turns nichrome wire surrounded 
fire brick plates reduce heat losses. 


the liquid agitated during growth, ripple mark frozen along the 
top edge the interface. agitating equal time intervals, the 
approximate rate growth the crystal may determined from subsequent 
measurements the position the corresponding ripple marks. The crystals 
were grown air, and the oxide which formed the liquid surface was 
skimmed off with flattened glass rod. 


During the growth number crystals rates controlled uniform 
movement the furnace, the furnace was suddenly pulled away and the end 
the boat dropped, allowing the remaining liquid run down the boat away 
from the solid metal. was found that slow rates growth (less than 
per minute) the interface was essentially plane and perpendicular 
the direction growth. higher rates, the interface became progressively 
more concave, owing the higher thermal conductivity the graphite. 
case, even rates 2.5 cm. per min., were dendrites observed. 


However, the furnace was completely pulled away from the boat after 
solidification had started, and the crystal allowed freeze rate governed 
only the rate which heat was lost the surrounding atmosphere, was 
found that large dendrites could obtained. This procedure was therefore 
adopted for this investigation, and all the crystals illustrated were grown 
this manner. The rate growth measured the top surface now appeared 
accelerate roughly exponentially until the first dendrite peak appeared 
the surface. Following this, the rate growth could not defined since the 
nominally metal consisted partly solid dendrites and partly 
liquid until solidification was complete. 


The seed crystals used were usually cm. long and were the full width 
the boat, ensure symmetrical heat flow the horizontal plane. 
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Bicrystals were grown essentially the same manner that the single 
crystals. the place the one wide seed crystal, two narrow seeds the 
desired orientations were placed the boat side side and separated 
short graphite insert. crystal was then grown uniform rate simultan- 
eously from both small seeds point well past the end the insert. The 
result was usually bicrystal which was the full width the boat. This was 
subsequently used the seed for crystals grown the usual manner. 


Measurement Temperature 


small thermocouple was inserted horizontally into the molten lead through 
holes made the side the boat. The temperature the metal, indicated 
the thermocouple, was then recorded the crystal grew. 


The thermocouple was made silver soldering copper and constantan wires 
having diameters 0.0048 in. The copper wire had glass fiber wrapping 
which satisfactorily insulated from the bare constantan wire for the tem- 
perature range encountered. The protecting tube surrounding the wires was 
made glass and had outside and inside diameter 0.9 mm. and 0.5 mm. 
respectively. 

Preliminary temperature measurements with Leeds and Northrup Type 
potentiometer showed that was necessary immerse the thermocouples 
least mm. into the liquid obtain consistent results the freezing point 
lead. With sufficient immersion, the measured temperatures are believed 
correct 0.1 degree Centigrade. 

After calibration, the thermocouple was placed 1.5 mm. above and parallel 
the bottom the graphite boat with the junction 2.5 cm. from one side wall. 
Molten metal was then added the boat and crystal mm. thick was grown 
the usual manner from seed the end the boat. Most the resultant 
e.m.f. the immersed couple during growth was balanced with the potentio- 
meter kept fixed position. The remainder was recorded Leeds and 
Northrup Speedomax recorder having full scale sensitivity mv., giving 
temperature measurements believed correct 0.2 degree Centigrade, and 
temperature differences taken from one curve, correct 0.1 degree Centigrade. 


series cooling curves was then obtained for one crystal which was re- 
grown number times. each case, the point which solidification started, 
with respect the fixed couple, was changed, effectively changing the position 
the couple the crystal. The necessary amount metal was added keep 
the total amount metal solidifying each case constant. 


attempt was made grow crystals thicker than mm. order de- 
crease the ratio thermocouple thickness specimen thickness; this would 
reduce any disturbance the system due the presence the thermocouple. 
However, the attempt was abandoned, since was found that stray crystals 
tended form the thicker crystals when they were reheated, even with very 
careful handling. The possibility introducing number thermocouples 
various fixed positions the crystal was also considered. was not adopted, 
since the presence number couples could appreciably change the normal 
thermal conditions during crystal growth. 
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Experimental Observations and Discussion 


Dendrite Structure and Manner Growth 


typical example the dendritic structure obtained shown Fig. 
which photograph decanted single crystal lead viewed from above 
and one side. The crystallographic orientation indicated the three 
arrows, each arrow being parallel <100> crystallographic axis cube 
edge, since lead has face-centered cubic structure. Growth proceeded from 
the left from point outside the photograph. 


The dendritic structure seen consist series parallel ‘‘dendrite 
each row corresponding the upper half the dendrite skeleton drawn 
Fig. Taking the primary stalk the skeleton being the horizontal 
plane the bottom the crystal, the vertical secondary and horizontal ter- 
tiary stalks one row are clearly seen, growing direction corresponding 
the diagram. The tertiary branches growing along the surface the liquid 
branch even further shown Fig. which photograph showing the 
details part the top surface vertical stalk. some cases, the larger 
vertical stalks have been disturbed the decanting process and are tilting over, 
the smaller ones remaining perpendicular the primary stalk. 


Fig. the secondary branches are shown progressively decreasing length 
and thickness one proceeds from the center the outer ends the primary 
stalk; this, however, not the case for the actual dendrites shown Fig. 
examining the vertical stalks, for example, apparent that relatively 
large number stalks begin grow, but only few reach the top surface where 
the vertical growth arrested. Nor the vertical stalks which develop 
grow equal intervals shown Fig. would also expected that, since 
the secondary stalks grow from the bottom the boat, the tertiary branches 
should larger and thicker the closer they are the bottom; this not the 
case here for the larger stalks. 


The manner which the dendrites grow indicated the three crystals 
shown Fig. These crystals were grown under similar conditions, but were 
decanted 50, 52, and sec., for crystals (a), (b), and (c) respectively, from 
the instant which the furnace was removed. The upper edge the 
liquid interface crystal (c) had advanced mm. further from the point 
which solidification started than that Comparing the three crystals, 
seen that the dendritic rows crystal (a) are shorter and thinner than those 
the other two. Their secondary branches are very small, showing further 
branching, and the spacing between rows smaller than that the other two. 
The photograph does not show some the smaller rows for this crystal. The 
same lesser degree may said the dendrite size and spacing with 
respect (c). The manner which the dendrites grow therefore appears 
agree generally with the description given Carpenter and Robertson, the 
primary stalk growing first, followed the secondary and tertiary stalks, 
some which progressively increase size. noted before, this not true 
for the vertical stalks which reach the top surface, since their growth arrested 
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Fic. decanted lead single crystal viewed obliquely from above. The arrows are parallel 
the crystallographic axes. 


Fic. the top branches vertical dendrite stalk. 
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Crystals decanted progressively later solidification. 
the dendrite pattern for two crystals grown under similar conditions. 
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Fics. AND 
The arrows are parallel their crystallographic axes. 


PLATE 


Fic. view decanted crystal orientated shown. 


Fic. Dendritic pattern bicrystal whose components are orientated similarly those 
Figs. and 


Fic. view crystal showing the transition from striations the top jagged 
boundary zones. 
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the surface and their branches grow much larger the top. the same 
time that the dendrites are growing, the general concave interface 
sweeps along covering the rows initially formed. 


Direction Growth 


The relation between crystallographic orientation and direction dendrite 
growth further shown Fig. which plan view crystal having the 
same orientation that Fig. The mutually perpendicular directions 
growth primary, secondary, and tertiary branches clearly seen, although 
two vertical stalks have fallen over the decanting process. The crystal shown 
Fig. was grown under similar conditions those above except that the 
direction now parallel the boat axis instead the [100] direction. 
The dendrite rows and smaller branches grow 45° the general direction 
growth, but still parallel the cube edge the crystal lattice. ensure that 
growth conditions were the same for both crystals, bicrystal was grown and 
the result shown Fig. The crystal the right has the same orientation 
that Fig. that the left the same that Fig. therefore appears 
that the direction dendrite growth depends only the crystallographic 
orientation the parent crystal and not the direction solidification. 


Crystals other orientations were grown and decanted, and the results con- 
firmed the above statement. Fig. shows end view crystal whose 
orientation such that one axis lies along the longitudinal crystal axis, the 
others being approximately 45° the horizontal plane. expected, the 
secondary stalks now grow correspondingly 45° the horizontal plane. 
For the cases which cube edge was not parallel the longitudinal crystal 
axis, for example the case for the [111] direction parallel this axis, 
dendrite rows were not observed, but stalks were observed near the completely 
solidified crystal, projecting directions corresponding the cube edges 
the crystal. should noted that crystal (c) Fig. has not quite the same 
orientation that and This accounts for the small difference 
dendritic row direction between and the other two. 


all cases observed, the dendrite direction appeared correspond within 
two degrees the crystallographic direction, this value being the limiting 
accuracy for the experimental procedure adopted. However, when crystals 
similar that Fig. were examined more closely, was found that all the 
dendrite rows were not exactly parallel, some rows diverging from their neigh- 
bors much three degrees. the same time, when the crystals were 
deeply etched 10% nitric acid solution, was found that all the crystals 
did not have the same orientation seen from the different optical reflective 
properties each section and verified back-reflection X-ray photographs. 
This phenomenon will discussed more fully the next section. 


determine the effect orientation differences direction dendrite 
growth, six bicrystals were grown which one component the bicrystal was 
related the other simple rotation about the vertical axis, Fig. 
Measurements were then made the position the dendrite rows with low 


q 
q 
q 
q 
q 
q 
q 
q 
q 
q 
4 
| 
q 
q 
4 
| 
q 
4q 
q 
q 
q 
q 
q 
3 
: 
dq 
. 


388 CANADIAN JOURNAL PHYSICS. VOL. 


powered microscope and travelling stage, and series X-ray photographs 
were taken traversing the crystal. The divergence the dendrite rows two 
zones, one each component the bicrystal, could then related the 
orientation difference the two zones. The results are listed Table 


TABLE 


Bicrystal Dendritic row divergence Orientation difference 


The accuracy measurement the direction dendrite row limited 
the round tips the vertical stalks which are lined the cross hairs 
the microscope. However, these experiments, this error less than the 
probable error the measurement the relative orientation the zones, 
estimated one degree. This due the size the Laue spot the film, 
the X-ray spot falling the specimen being 0.7 mm. diameter. 


may, therefore, concluded that, within experimental error, the direction 
dendritic growth uniquely defined the crystallographic orientation 
the and independent the direction heat flow. The presence 
zones slightly different orientations could account for the fact previously 
noted that the dendrite rows single crystal are not quite parallel. 


Crystal Substructures High Rates Growth 


slow rates growth substructures have been described for 
tin single crystals Teghtsoonian and Chalmers (6) who have also referred 
the existence similar structures lead. They report that striations ap- 
pearing the top surface crystal consist series bands parallel 
the longitudinal axis the crystal, with average width about milli- 
meter. The crystallographic orientation the striations with respect any 
one reference striation defined small simple rotation about the longi- 
tudinal axis the crystal. higher rates growth, the average width the 
striations decreases and they converge, since they grow perpendicular the 
interface which becomes concave for crystals grown the present 
manner. They can seen Fig. the striations appearing the completely 
solidified portion the crystal the left-hand side the photograph. 


When crystals lead showing dendritic growth were deeply etched and 
examined, striations were observed the bottom the slowly grown seed 
well the top. However, the rate solidification increased, the appearance 
and nature the striations changed shown Fig. This the bottom 
view etched crystal similar that Fig. which some the stria- 
tions the seed crystal appear the top the photograph. can seen 
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that few striations widen considerably the expense the others. Further, 
the boundaries between them are longer straight lines but become jagged 
and not remain parallel the average the longitudinal axis the crystal. 
X-ray examination these jagged boundary zones shows that each zone 
uniform orientation, but neighboring zones appear differ relative 
orientation small amounts, three degrees. There appears 
common axis rotation relating the zones any one crystal. The average 
direction the jagged boundaries tends parallel the dendrite rows. 


Looking the top surface the crystal, jagged irregular boundaries corres- 
ponding those observed the bottom are seen the decanted portion 
the crystal. This shown Figs. and 11. Fig. the plan view the 
top surface crystal illuminated bring out some the zones. Four zones 
can seen the lower half the crystal. Fig. photograph the 
bottom surface the same crystal, showing the same zones. shown 
Fig. 10, the jagged boundaries run between dendrite rows appearing 
compromise between the horizontal side arms projecting primary 
stalks the rows either side the boundary. all the specimens 
examined, dendrite row ever crossed one these boundaries, although 
some cases rows terminated when they encountered one. Also some cases 
where large dendrites grew from the lip from the side the crystal toward the 
center, before the dendrite rows had grown that far, the boundary turned 
sharply and followed the direction these lateral dendrites. The boundaries 
the decanted portion could clearly distinguished very nearly the 
upper edge the concave interface the crystal where the regular 
striations appear. one case, the crystal was electrolytically etched weak 
solution nitric acid. This resulted the striations’ receding approximately 
mm. toward the seed after 0.27 mm. material had been removed from the 
top surface, bringing out zones which were continuations those observed 
the decanted portion. Further etching did not materially affect the end position 
the striations. 

Jagged boundaries were also observed the bottom surface bicrystals. 
The degree irregularity does not appear affected the difference 
orientation the crystals each side the boundary, but appears more 
irregular the the crystal that solidify last where the rate growth 
probably greatest. 


tentative explanation the development these jagged boundary zones 
advanced terms dendritic growth. rates growth greater than about 
one centimeter per minute lip observed the decanted surface, projecting 
above the bottom the boat the foot the nearly plane vertical solid— 
liquid interface, indicating that the part the liquid adjacent the graphite 
boat slightly cooler than the remainder. Therefore, possible that very 
small flat ‘‘fern dendrites could grow along the bottom the boat taking 
the orientation the striations from which they originate. However, 
result the slight differences thermal contact between the liquid and the 
graphite boat, some dendrites may grow ahead others, the former being the 
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basis the zones. Once this condition has been reached, the lateral branches 
the leading dendrites will grow front the slower dendrites preventing 
their further growth. also possible this stage that convection currents 
the liquid adjacent the solid, vibration, could bend the fine dendrites, 
accounting for the small orientation differences neighboring zones. When 
the lateral branches the remaining dendrites meet, one would expect 
jagged boundary compromise between the two, observed, depending 
only the outline the dendrites and not their relative orientations. 


noted before, striations grow perpendicularly the inter- 
face; therefore, would expected that those growing along the bottom 
the crystal would tend swing towards the top order that they remain 
perpendicular the lip surface. new striations did not form, the metal below 
the striations would then take the orientation the solid defined the 
small dendrites. Therefore, since the lip becomes progressively larger growth 
accelerates, the zones forming from the small dendrites along the bottom would 
expected become thicker until the time the decanted portion the 
crystal reached, they would extend through the entire thickness the 
crystal. Since the very top the interface remains essentially 
vertical, shown Fig. the striations would expected extend much 
farther along the crystal the top than the bottom. They should disappear 
near the decanted portion when layer the top surface removed etching, 


was observed. 

The fern type dendrites referred have been observed number cases 
growing from the solid along the top the liquid exposed the air. the 
name implies, they are flat with closely spaced side branches and often project 
several centimeters ahead the interface into the liquid. example shown 
Fig. which the fern type dendrite has fallen the bottom when the 


liquid under was poured away. 

The effect thermal contact between the liquid and the graphite boat 
rate growth shown markedly Fig. 13. Normally, when metal added 
the boat and melted, layer oxide forms both the top surface the 
liquid and the bottom surface between the liquid and graphite boat. The 
top layer oxide removed, but the bottom layer usually remains, resulting 
poorer thermal contact between liquid and boat. the crystal Fig. 13, 
the oxide was removed along strip running down the center the crystal, 
rubbing the boat surface under the liquid with flattened glass rod. 
seen that the rate solidification above the strip was greater than that the 
remainder, shown the convex portion the general concave solid—liquid 
interface. removal the oxide affects the rate solidification the extent 
shown, then seems feasible that slight differences temperature could exist 
along the boat surface produce the advance several tiny dendrites ahead 


their neighbors. 
Northcott (3) notes his examination dendritic structure brought out 
polishing and etching section copper alloy casting that irregularities 
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the structure were present what appeared single crystal. X-Ray 
examination showed small differences orientation between these regions. 
appears that this irregularity could due the presence zones similar 
those described above. 


Supercooling 


Some the cooling curves obtained the manner previously described are 
shown Fig. 14. Their corresponding position the crystal given the 
distance between the thermocouple and the point which solidification 
started, Table II. This table includes the supercooling indicated the 


TEMPERATURE 


10 
TIME 


Fic. 14. Cooling curves for various positions the crystal. 


TABLE 
Curve Supercooling, Dendrite appearance 
6.5 Well developed—as the decanted 
5.5 1.2 portion Fig. 
4.5 
3.5 0.5 Partly developed—as (0) 
0.0 None 
0.5 0.0 None 
| 


curves and the dendritic pattern observed the region surrounding the thermo- 
couple for crystals grown under similar conditions. 


From the last two columns Table II, appears that supercooling 
necessary condition for dendritic growth high purity lead. also appears 
that the development the dendrite branches related the degree super- 
since for thermocouple positions giving the same supercooling (Curve 
the dendrite development generally the same; for less supercooling 
(Curves and the development correspondingly smaller. 


has previously been noted that dendrites were not observed when crystals 
were grown uniform rates. would therefore expected that under these 
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growth conditions the liquid ahead the interface should not 
supercooled. This was found so, within the sensitivity the measuring 
apparatus. 

Summary 


For high purity lead, concluded that dendrite skeleton forms from 
small single stalk projecting into the liquid ahead the general 
interface. The secondary branches the skeleton grow from and perpendicular 
the primary stalk grows larger and thicker. Some these secondary 
branches become larger growth proceeds and develop tertiary branches 
which, turn, branch again and on. 


All the stalks grow parallel <100> directions defined the crystal- 
lographic orientation the dendrite. Their directions are not influenced 
the temperature gradients present the liquid which they grow. 


high rates growth, zones similar crystallographic orientation are 
observed the crystal, each zone having small difference orientation with 
respect its neighbor. Because this orientation difference, all the dendrite 
rows single crystal are not exactly parallel, since their direction defined 
the orientation the zone which they grow. 


believed that the orientation the zones determined the initial 
tiny dendrites which develop during growth and which some cases have been 
disturbed convection currents and vibrations the liquid surrounding 
them. The zones are characterized jagged irregular boundaries which appear 
result compromise between the side branches the two dendrites 
growing either side the boundary. They are not observed crystals grown 
slow rates growth. 


From temperature measurements, appears that supercooling necessary 
condition for the dendritic growth high purity lead, and that the size and 
spacing the dendrite stalks function the degree supercooling the 
liquid which they form. 
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THEORY MULTIPOLE 


WALLACE 


Abstract 


develop systematic analysis the radiation from given oscillating 
system charges and currents, without any approximations. Using simple 
vector identity, the vector potential separated into irrotational and solenoidal 
parts. The field may expressed terms the latter alone. similar vector 
identity involving the operator grad (the rotation operator) permits the 
separation the field into parts which the radial components the electric 
and magnetic field, respectively, vanish. The energy flux, energy density, and 
angular momentum density may each case expressed terms the angular 
operators Expansion the eigenfunctions these operators, the spherical 
harmonics, corresponds the separation into electric and magnetic multipoles 
all orders. Introduction spherical enables exhibit 
these radiations terms natural multipoles (derivatives All calcu- 
lations are carried out without restriction size radiating system relative 
wave length, the induction well the radiation region. 


Introduction 


There exists extensive literature the theory multipole radiation 
(2-9). Some this concerned with specifying the form the radiation 
fields for the various radiations, without reference the source the radiation. 
The most elegant treatment undoubtedly that Franz Those papers 
that deal with the radiation from given system oscillating charges and 
currents are general not free approximations based the assumption 
wave length large compared with the dimensions the radiating system, 
are restricted the distant field only, the lower multipole orders.* 
The present paper represents attempt give general treatment, using 
only elementary mathematical techniques. The aim proceed far 
possible with the analysis the field merely into and 
parts, and express the energy, angular momentum, etc., each. Expan- 
sion these quantities spherical harmonics then corresponds immediately 
separation into multipole orders. 


Representation the Field Vector Potential 
shall suppose that the source the radiation consists charge dis- 
tribution p(x’, 2’, po(x’, y’, and current distribution J(x’, y’, 


received April 1951. 

Contribution from the Department Mathematics, McGill University, Montreal, Que. 

However, Dr. Blatt the University Illinois has been kind enough show the author 
alternative general treatment, the spirit the paper Franz, which will appear forth- 
coming text book (on nuclear physics) Weisskopf and Blatt. 

the time variation not simple harmonic, the charges and currents must analyzed 
into Fourier components; the field then superposition oscillatory components with distri- 
bution frequencies. 
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becomes this case 


terms and the scalar and vector potentials the field may written 


where being the point which the potentials are being evalu- 
ated and the variable integration ranging over the distribution. With 
the given distributions (2.2) may written 


where w/c being the wave length the radiation. 


Consider now the vector identity 


where the gradient operator with respect the primed co- 


ordinates. Letting this and using the fact that 
get for 
ikR 
may now use the fact that 
ikR ikR 
(2.6) 


and integrate the first term parts obtain 


ikR 


The two parts are respectively irrotational and solenoidal. The first 


therefore contributes nothing the magnetic contribution the 


electric field, may written terms the charge density 


virtue the continuity relation (2.1); dV’ and therefore cancels 


Thus, the whole field may derived from the vector potential 


If, addition currents, there oscillating density magnetization one must add 
mechanical problems. 
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Separation into and Radiations 


The key the separation (2.8) into and radiations 
the introduction the operator” 


(3.1) 


This operator the same the quantum-mechanical angular momentum 
operator within (complex) factor. may also envisaged follows: 
the point rotated with angular velocity about the origin. 
follows that expressed polar coordinates and expanded spherical 
harmonics, the radial dependence unaffected the operator i.e. 
operator the angular function only. useful result, which follows immedi- 
ately, the following: function the angle between the position 
vectors and r’, then 


(3.2) 


This follows from the fact that rotation about any axis induces the 
same change and therefore equal and opposite rotation 


The splitting the vector potential (2.8) now follows the use vector 
identity similar (2.4), but with L-operators replacing the gradients: 


L.L'= 
and may write (3.3) symbolically 


where the inverse operator L?. Let multiply the right 
Now introduce function such that 


eikR eikR 


where designates the average over-all directions. 
AV 


using the known expansion for spherical harmonics, find that 


being the Bessel function” order and (kr) the spher- 
ical Hankel function the first kind (diverging wave). 
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Now since 
ikR 


get, letting (3.4) operate and integrating, 
ikR 
Therefore 


the second, 
LV’, 


the second becomes 


c 


This, shall see, the potential the magnetic multipole radiations, while 
gives that the electric. 


shall now transform more convenient form. Equation (3.8), 
for the last may replaced Then may write 


etc. Using the commutation rules for the components 
and the fact that this becomes 


may interchange the order the and operators the second term 
using the commutation relations the form: 


find then that 
(J.N)Lz = L.(J.N) + (J X N); 


may noticed that the operator the square bracket (3.4) automatically eliminates 
spherically symmetrical components the operand. Since not necessarily zero, 
conceivably the operator could restore terms this type. However, the omission terms 
what follows (e.g. the function direct consequence the earlier removal the irro- 
tational part the vector potential. 
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from which follows that 
the remaining terms cancel since 


still more symmetric form may obtained use the fact that 
that may write finally 


The identification (3.10) and (3.14) the potentials for the magnetic 
and electric radiations respectively may made writing down the cor- 
responding fields: 


Cc c 


c c 


But now that the radial components and are zero. 
the radial magnetic field zero, describe the radiation electric 
the radial electric field zero, magnetic multipole. 

The integrals occurring (3.10), (3.14), and (3.15) may written 
more suggestive form. For 


Using this and (2.1) get 


The second term bears the first ratio a), where the radius the 
radiating system, and consequently negligible cases which the dimen- 
sions the radiating system are small compared with the wave length the 
radiation. 


*See Appendix. 
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From (3.15) the fields the case electric radiation are 
and the magnetic 


The energy flux given 
P = Re + x (H,* + 
us 


Where denotes “real part of’’, and designates complex conjugate. 
But may verified directly that 
div (Ey x H.*) = — div xX H;*), (3.20) 


and the contributions the energy from the electric and magnetic radiations 
are independent each other. 


Radiation 


The energy flux electric radiation over sphere radius 
us 


where the integral over solid angle. Now from (3.17) 


This enables write (4.1) 
us 


where designates part The integrand (4.2) gives the 
radiation pattern, i.e. the radial energy flow function position general. 


evident from this that the various electric multipole radiations are 
independent, since the spherical harmonics terms which shall expand 
are eigenfunctions Thus the respective contributions the flux are 
independent. 


Taking the imaginary part rather than the real part (4.2) could find 


(1) the difference between the mean electric and magnetic field energies. Since 
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The angular momentum the radiation also easily found; the angular 
momentum density 


Thus the expansion spherical harmonics will also expansion the 
eigenvalues angular momentum. Furthermore, the angular momentum 


density goes asymptotically that there the same angular momentum 
any shell dr). Therefore, the angular momentum effectively all 


the radiation field.* 


The Radiation 


The results here parallel very closely those the preceding section. The 
energy flux may written 


us 


Since operates only the radial part G2, expansion multipoles again 
results complete separation the terms various orders. Also the 
angular momentum density comes out 


Multipole Expansion 


Let now introduce the spherical harmonics” (10) 


and make use the fact that 


‘The indicates dummy indices summation and the angle that 
between the directions and may now use the expansion (3.6) for 


*However, there would angular momentum the fields were both without radial com- 
ponents. The decrease the radial component with compensated the increase the volume 
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substituting from (6.2) for (cos obtain expansions harmonics 
(corresponding the various multipole orders*) for and the fields then 


follow from (3.18) and (3.19). 


The radiation intensity may therefore written 


(6.6) 
where 

and the refer different sets dummy indices. 


where the sum the different products Kronecker 6’s involving 
one index from and one from and the involve 
least two indices from the set (/) These latter terms give con- 
tribution since, their definition, with repeated index vanishes. 
Thus the flux radius 


where the scalar product and its complex conjugate. For 
large the radial factor approaches unity. 


Similar calculations with the terms lead 


follows that the flux the magnetic radiation has intensity 


2 


should noted that, when the expansion carried out this corresponds 
dipole radiation (electric magnetic), quadrupole, etc. 


. 
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where the last term again becomes unity large distances. must em- 
phasized that the formulae (6.9) and (6.12) are valid without any approxi- 
mations concerning the relative magnitude the dimensions the radiating 


system and the wave length. 

are interested the components angular momentum, the tensor 
approach probably longer the natural one. Rather, seems desirable, 
one calculate the z-component angular momentum, use spherical 
harmonics terms polar coordinates with the z-axis polar axis; i.e. 
expand terms the eigenfunctions and L,. 


also some interest calculate the energy density. For the 
type radiation this 


Integrating this over volume, integrating parts and using the fact that 


and for the radiation 
2 
av. (6.14) 
now expand terms “conventional’’ spherical harmonics 


designating the coefficients for brevity, evident that the con- 
tributions the z-component angular momentum are 


while the corresponding contributions the energy are 
Thus have the general result 
7( lm) 


For completeness, list the expansion coefficients 
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Both the total angular momentum and total energy will course finite 
only the system has been radiating only for finite period time. 


Appendix 


may write (cf. Equation 3.5) 
e kR 


where the inverse operator L”. Then 
ikR 


1-2 
L’ 


now necessary only know the commutation rule between and L”: 

and that between and 

Multiplying the left, get 

Putting this (A.2) and operating with another gives the required result: 


(A.3) 
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VARIATION METEOR ECHO RATES WITH RADAR 
SYSTEM PARAMETERS! 


Abstract 


Observations made with crossed-polarization radar system systems not 
support the suggestion that the ionized meteor trail may act strong 
polarizer the incident radio wave. Experiments have been carried out 
determine the variation normal meteor echo rates with transmitter power, 
antenna gain, and radio wave length, and all confirm Lovell’s scattering formula, 
provided that account taken the effective broadening the scattering 
pattern the meteor trail with increasing wave length. The limiting sensitivity 
the 9.22 200 kw. radar determined about 9th magnitude. During 
strong visual shower the observed increase visual rates and low-power radar 
rates, compared high-power radar rates, explained assuming that the 
magnitude distribution the shower meteors differs from the normal nonshower 
distribution. 


Introduction 


From earlier work Blackett and Lovell (2), Lovell and Clegg (8) have 
developed theoretical relation between the electron density meteor trail 
and the radar system parameters, assuming the echo signal produced 
scattering from long line free electrons. When the numerical values for 
the physical constants are inserted this formula becomes, 


2,3 


3 \1/2 
5.5 (=) electrons per cm. 


received echo power watts 

number electrons produced per cm. path 
antenna gain 

distance meteor trail cm. 

peak transmitted power watts 

wave length transmitter cm. 


Lovell (8) demonstrated that the predicted variation with was agree- 
ment with the observed variation over the range 1.4 8.3 m., for indi- 
vidual echoes the type that could generally associated with visual 
meteors. earlier paper (7) showed that the number observed echoes 
falling given interval about selected value was approximately inver- 
sely proportional From theoretical considerations, Herlofson (5) indi- 
cated that Ist magnitude visual meteor moving km. per sec. may create 
electrons per cm. path and 6th magnitude meteor about electrons. 
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However, Lovell (9) has suggested that the formula may not applicable 
the brighter meteors detected longer wave lengths. The question first 
arose comparing the rates meteor echoes 8.3 (and 9.2 
(11)) the rates echoes 4.2 and shorter wave lengths. There appeared 
sudden increase the hourly counts for wave lengths greater than 
about m., and Lovell (10) suggested that other scattering processes may also 
occur the longer wave lengths. interpreted the phenomenological theory 
McKinley and Millman (11), based their observations 9.2 m., 
implying that multiplicity echoes can arise from single meteor. Herlofson 
(4) recently proposed theory plasma oscillations the column ionized 
gas formed the wake the meteor. One consequence this theory was that 
under certain conditions obtaining the longer wave lengths the echo strengths 
could many magnitudes greater than would expected from the formula, 
Equation 


Alternatively, one might regard the ionized cloud left the meteor 
miniature ionosphere. the cloud assumed have effective echoing area, 
the received signal power may found from the conventional radar 
equation (14), 


€ (2) 


The volume density electrons must exceed the critical density for the wave 
length employed, otherwise echo will received regardless the size 
the echoing area. The critical density inversely proportional the square 
the wave length. the wave length increased fainter meteors should 
detected, both for this reason and also because the factor Equation (2). 
effect, the minimum detectable echo power should depend for this 
total reflection case, compared with for the case involving scattering from 
long line source free electrons. 


the course general program meteor observing, which has been 
progress Ottawa since the summer 1947 under the combined auspices 
the Dominion Observatory and the National Research Council, some experi- 
ments have been conducted 9.2, 5.4, and 2.8 m., with various transmitter 
powers and antenna systems. These have led the broad conclusions that 
Lovell’s formula applies reasonably well the majority meteors detected 
wave lengths long 9.2 However two new factors are brought for- 
ward which have bearing the variation meteor echo rates. One 
further decrease the number detectable echoes the wave length 
decreased this effect consequence the geometry the diffraction pattern 
associated with the scattered echo from the meteor trail, and has not been 
included Equations la, 1b. The other distinct difference the size 
distribution meteors belonging strong shower compared with the 
distribution the normal background meteors. The following sections, 
order, describe experiments designed test Herlofson’s plasma theory, the 
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variation echo rates with antenna gain, the variation rates with trans- 
mitter power, the variation rates with wave length, and the general con- 
clusions. 

Polarization Effects 


Herlofson’s theory (4) plasma oscillations suggests that, the volume 
density electrons the meteor trail sufficient cause the dielectric 
constant the medium become negative, anomalous effect should 
observed which dependent the polarization the incident electromagnetic 
wave relative the axis the meteor trail. the electric vector parallel 
the trail the usual scattering processes apply, defined Equation (1). the 
electric vector normal the ionized column, plasma resonance oscillations 
are presumed occur, with the result that the scattering should become several 
powers stronger than for the parallel case. Herlofson adds that, owing 
this resonance phenomenon, ionized column would act combined filter 
and polarizer for the scattered radiation. would return with great strength 
radiation with such low frequency that the dielectric constant negative 
the center, but would respond feebly radiation higher frequencies. 
the favorable frequencies only the wave component with electric vector normal 
the cylinder would scattered with appreciable intensity. 


the incident electric field vector angle the meteor column 
may resolved into components cos and sin which are respectively parallel 
and normal the cylinder, and shall assume from the above discussion 
that only the normal component will contribute significantly the received 
echo signal. Our experience (11) has shown that the great majority meteor 
echoes observed 9.2 occur near the point where the meteor path 
normal the line sight: the brighter meteors can certainly detected over 
much wider aspects but their number statistically insignificant. For the 
majority echoes, therefore, the angle will lie the plane the incident 
wave front, which simplifies the geometry the situation. observational 
periods are selected when major meteor showers are active reasonable 
assume that the sporadic meteor paths will random directions; that is, 
there will fairly random distribution the angle 


Case two similar but independent radar systems, one operating 
with its transmitting and receiving dipoles oriented and the 
other Each antenna half-wave horizontal dipole mounted 
quarter-wave above ground screen reflector. The wave lengths differ slightly 
avoid mutual interference. The electric vector the N/S transmitter, ex- 
pressed volts per cm. the target distance, will produce after reflection 
from the meteor column scattered field strength proportional sin the 
vector being normal the column. This vector will have component pro- 
portional sin the direction the receiving dipole the N/S system 
that the resultant voltage produced the input the receiver will pro- 
portional sin? Similarly the receiving dipole the E/W system will 
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excited field proportional cos? the receiver outputs vary linearly 
with field strength the output signals will respectively proportional 
and cos? and the ratio echo voltages the two systems will propor- 
tional tan? 


Case one radar transmitter with its dipole oriented N/S, and 
two receivers both tuned the transmitter frequency, one having its dipole 
N/S, and the other E/W. The N/S receiver will produce echoes proportional 
sin? before, but the E/W receiver output will now proportional 
sin cos and the echo voltage ratios will proportional tan 


One the basic requirements the two cases described above that the 
electric vectors the N/S and the E/W systems shall right angles 
each other the meteor target, and this condition holds accurately only the 


RANGE CLASS 
KM. 100-120 KM. 


HORIZON 


Fic. projection the sky. The shaded areas show the meteor zones corresponding 
the selected range classes. The curves, marked 2.5° 30°, indicate the departure from ortho- 
gonality the N/S and E/W systems. 
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zenith and along the two great circles passing through the zenith and projected 
axes the dipoles. 


Fig. zenithal equidistant projection the sky which have been 
drawn curves indicate the departure from orthogonality the electric 
vectors the two systems, the selected values being 2.5°, 5°, 10°, 20°, 30°. 
assumption that the part the M-region lying between km. and 100 km. 
height above sea level more productive meteor echoes (11, 13), heavily 
shaded areas have been drawn show the zones the sky from which most 
the echoes given radar range class may expected occur. The lightly 
shaded areas denote the height regions, 80-90 km., and 100-110 km., from 
which fewer echoes are expected. Because the records the radar 
systems carry km. calibration lines convenient consider range classes 
steps km. The sky coverage for the range class 80-100 km. shown 
the upper left quadrant only Fig. that for the range class 100-120 km. 
the upper right quadrant, and on. will immediately apparent that 
the range class 80-100 km. the best one consider for our purposes the 
departure from orthogonality for most the echoes should less than 2.5° 
and for very few great 10°: The range class 100-120 km. embraces sub- 
stantial region where the departure excess but the statistics for this 
class should still some value. For the range classes 120-140 km. and greater 
the simple analysis breaks down, and, fact, for the long range classes low 
angles the polarizations the two systems are not far from being parallel for 
great percentage the echoes. 


further reason, and very practical one, for rejecting data from the long 
range classes was the marked difference the antenna coverage the N/S 
and the E/W systems low elevation angles. There significant difference 
the theoretical radiation patterns dipole backed reflecting screen, 
depending whether the measurement made the plane the dipole 
normal it, though this difference importance only low angles. Further- 
more, under the actual experimental conditions, several towers, buildings, 
poles, and guy-wires the neighborhood the transmitting and receiving 
dipoles caused additional distortions the antenna patterns low angles. 
The distortions were not the same for the two systems because the different 
effects the polarizations the various reflecting objects the immediate 
vicinity, particularly the antennae were separted few yards. un- 
likely, though, that any serious difference coverage existed high angles. 


For study Case two transmitters were used which were similar but not 
quite identical. The E/W transmitter used two tubes the oscillator with 
power output about 200 kw. peak, whereas the N/S oscillator used four 
tubes and the measured relative power output was 40% greater. Both trans- 
mitters were keyed parallel from the same modulator with the same pulse 
width, pulse voltage, and recurrence frequency. avoid mutual interference 
the N/S system operated 9.22 and the E/W system 9.52 With 
this separation each receiver had rejection ratio more than db. for un- 
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desired desired echo signals. comparing echo strength the two systems 
small correction factor should applied. Employing Lovell’s formula 
(Equation the difference the systems amounted power ratio 
1.25 db. favor the N/S system, which was negligible comparison 
with the probable error db. involved estimating most the ratios 
echo amplitudes the two systems. Equation may not the correct 
relation use, priori, but the differentials power and wave length are small 
and tend cancel, hence unlikely that any other reasonable formula 
(Equation (2) for example) would show any significant difference. 


Two identical receivers were used, and the outputs displayed two cathode 
ray tubes mounted side side and photographed simultaneously camera 
which the film was drawn vertical plane continuously past the open lens 
in. per minute. The same horizontal range sweep and calibration marks were 
applied each C.R. tube and the receiver outputs were fed the intensity 
grids the C.R. tubes. During the runs the N/S and E/W antennas were 
interchanged frequently between the receiver-display systems eliminate 
any discrepancies performance, though care was taken equalize them 
initially far possible. The video output each receiver passed through 
adjustable clipper circuit which could set limit the amplitude strong 
signals any given level, before application the C.R. tube. 


these experiments, the background noise level output the receivers 
was set the optimum value for operation the displays, and the clipping 
level was set twice the mean noise level, measured volts A-scan 
monitor C.R. tube. The receiver band widths were greater than the minimum 
required for the transmitter pulse width (20 that weak echoes were 
easily detectable virtue the shape the echo spot compared the 
finer-grained noise spots. the average, echoes half the mean noise voltage 
was clearly resolvable the film. Thus, the working range the receiver 
between the minimum detectable signal and signal just reaching the clipping 
level was voltage, db. Signals above the clipping level produced 
black shadow, increasing range the display proportion the signal 
strength, due the video overshoot effect the clipping circuits. The vari- 
ation this shadow length with input signal voltage could calibrated fairly 
accurately the point where the receiver began saturate prior the 
clipping circuits. arbitrarily defined db. the absolute power level the 
weakest detectable signal; therefore db. was the absolute power echo 
just reaching the clipping level. was found that db. represented the 
strongest signal that could measured before the receiver characteristic be- 
came nonlinear. Signals below db. were not detected, course, and signals 
above db. could not accurately measured, but there were very few above 
the db. level. 


For statistical purposes was therefore convenient divide the working 
range receiver into two parts, the first extending from db. and 
including all signals from the weakest perceptible signal one just touching 


McKINLEY: METEOR ECHO RATES 409 


the clipping level, and the second including all echoes over the db. absolute 
power level. These signal strength groups may called absolute classes, based 
the arbitrary zero defined above, distinguish them from the classes 
which shall use classify the observational data. That is, shall endeavor 
list all the echoes falling either into relative class under db., wherein 
the difference the echo strengths the two displays less than db., 
into second relative class over db., which includes all the remainder. 


the reduction the records the echoes were considered pairs, using 
range the nearest km., time 1/10th second, and echo shape cor- 
relate the pairs. echo pair appeared both the N/S and E/W records, 
each with signal strength between and db., the pair was assigned the 
under class immediately. one echo was below the clipping level and the other 
above more careful estimate the relative signal strengths was required 
before classifying the pair either under over. both echoes were above 
db. the pair was classed under: few these should have been labeled 
over, where one both echoes the pair reached saturation, but the number 
this kind was statistically insignificant. 


When one member pair was missing the pair was classed over. this 
case, the observed echo lay between and db., the missing echo might 
have been just below detectability but still within db. the observed echo. 
The probable number such pairs which should have been marked under was 
means negligible, most the observed echoes lay close the db. 
absolute level, but was difficult determine the true proportion without 
making assumptions advance regarding the distribution echoes with signal 
strength. Without applying such correction should noted that this 
method analysis will put too many echoes the over class. Table sum- 
marizes the echo classification techniques. 


Case the theoretical voltage ratio signals the two systems 
expressed For db. difference echo power (4:1 field 
strength echo voltage), can lie between the value and 
Thus, the fraction echoes expected the under class the ratio 


TABLE 
CLASSIFICATION 


Absolute class observed echo Relative class echo pair 
Either receiver The other receiver difference signal strength) 

db. echo Over 

0-12 db. 0-12 db. Under 

Over db. echo Over 

Over db. 0-12 db. Under over (see text) 

Over db. Over db. Under 
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Expressed per cent, the simple filter—polarizer 
theory predicts that 41% should listed under and 59% over. Case 
where the voltage ratio proportional tan the corresponding angular 
ratio (tan~'4 whence follows that 69% should under 
and 31% over. now the other extreme and postulate that polar- 
ization effects whatsoever are introduced the mechanics reflection from 
the meteor trail, i.e., that the plane polarization the reflected wave re- 
mains the original plane, then follows immediately that for Case 100% 
the echoes should listed under and none over, while for Case none 
should under and 100% should over. another alternative assume 
that every echo circularly polarized then both Case and Case 100% 
the echoes should under. The and the polar- 
hypotheses are academic interest but may useful suggesting 
limits. 


The two independent systems (Case were operated 
from 1100 E.S.T. 1500 E.S.T. Nov. 26, 1948, and from 1700 E.S.T. Dec. 
1200 E.S.T. Dec. 1948. part the films was discarded because 
interference from local and distant sources made difficult times estimate 
relative signal strengths accurately. About seventeen hours satisfactory 
records were analyzed. The echo coincidences were noted pairs, according 
the procedure Table and the numbers under and over echo pairs are 
shown Table II, Row with the percentages given Row the over 
classes the numbers Type pairs are indicated brackets. Row contains 
the predicted percentages based the filter—polarizer theory, Row the per- 
centages assuming polarization effects, and Row the percentages for com- 
plete circular polarization. 


The observed percentages echoes the over class should reduced appreci- 
ably. Most the echoes with both members the pair showing are the under 
class, hence reasonable expect that large fraction the Type echoes 
should also the under class. However, even without attempting adjust 
the figures this respect, seen that the wide divergence between the 


TABLE 
ECHO DISTRIBUTIONS FOR CASE 


80-100 km. range class 100-120 km. range class 
Under Over Under Over 
Observed numbers 402 (44) 1837 322 (208) 
Uncorrected percentages 88% 12% 85% 15% 
theory 41% 59% 41% 59% 
Circular polarization 100% 100% 
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observed percentages over and under classes and the percentages predicted 
the theory indicates that the meteor trail does not act 
combined filter and polarizer the incident radio wave, least, not the 
extent that the component normal the trail several orders magnitude 
stronger than the component parallel the trail. would inadvisable 
conclude that these observations refute the possibility plasma oscillations 
occurring the meteoric cloud; they merely imply that the deductions made 
from the plasma hypothesis are not correct for the majority meteor echoes. 
the other hand the data presented far cannot said favor either 
the extreme hypotheses, since both predict 100% under Case More infor- 
mation required, some which may available from study Case II. 


The N/S transmitter only was operated with both the N/S and E/W re- 
ceivers tuned the transmitter wave length 9.22 (Case for total 
hr. the nights Mar. and Mar. 30, 1949. From these records six 
hours were selected which were satisfactorily free interference. Table III 
summarizes the observations, together with the various theoretical predictions 


for this case. 
TABLE III 


ECHO DISTRIBUTIONS FOR CASE 


80-100 km. range class 100-120 km. range class 
Over Under Over 
Observed numbers (26) 228 425 
Uncorrected percentages 40% 35% 65% 
Filter—polarizer theory 69% 31% 69% 31% 
Circular polarization 100% 100% 


Here, again the Type pairs, shown brackets, contain some that should 
have been listed the under class, though the proportion would less than 
the preceding case because the proportion known under class echoes less. 
correction were applied would tend narrow the gap between the 
observed percentages and the percentages predicted the filter—polarizer 
theory, rather than widen Case One point that the tabulated 
data does not emphasize that the majority the echo pairs both the 
under and the over classes were stronger the N/S receiver. fact, there were 
only three echoes six hours, both the 80-100 km. and the km. 
range classes, that appeared the E/W receiver and not the N/S receiver. 


‘The results Case should therefore regarded inconclusive far the 


theory concerned. 


Considering the data for Case and Case together appears evident 
that the meteor trail does not act with efficiency 
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several orders ten suggested Herlofson, though the observations 
not exclude the possibility that the effect might present but much smaller 
magnitude, e.g., factor five, that might applicable only 
few meteors under special circumstances. However, comparing the obser- 
vations with the polarization” percentages Table and Table III, 
seems probable that the returning wave largely polarized the incident 
plane (certainly not circularly polarized), with appreciable content 
random polarization. 


Variation Rates with Antenna Gain 


general rule the transmitter and the receiver the 9.2 radar system 
Ottawa have each been coupled separate single half-wave folded dipole, 
mounted quarter-wave above reflecting ground screen (rather than using 
the same antenna transmit and receive). For some observations Mar. 
and 10, 1948, three parasitic director elements were added each the 
antennae increase the directivity vertically upwards. was not practical 
measure directly the relative power gains the two 9.2 antenna systems, 
hence scaled-down models the antennae were built for wave length 1.5 
The antenna with three directors was determined have power gain 2.2 
over the antenna without the directors, measured point several wave 
lengths distant and vertically above the antenna. 


The original intention the experiments with the directive antenna system 
was look for meteor echoes that might occur well above the M-region (11), 
possibly the upper parts the E-region the F-regions. conclusions 
could reached because any evidence echoes that might have occurred 
heights km. overhead was masked the number echoes detected 
low angles the side lobes the antenna pattern slant ranges 
km. More recent experience has indicated that the existence appreciable 
number meteor echoes the higher altitudes very unlikely How- 
ever, the experimental data obtained with the high gain antenna may 
some interest here correlating echo rates with antenna gain. 


Mean hourly echo rates for each the 20-km. range classes from 320 
km. were obtained from representative six-hour period observing with the 
high-gain antenna, namely E.S.T. Mar. 10, 1948 (see Fig. dotted 
curve). The minimum 150 km. corresponds the first null the antenna 
pattern and the smaller maxima the side lobes. For comparison purposes 
similar data obtained with the low-gain antenna are shown plotted the solid 
line Fig. averaging the periods E.S.T. Feb. 1948, and 
0900 E.S.T. Apr. 22, 1948. Observations with the low-gain antenna should 
preferably have been carried simultaneously with the high-gain experi- 
ments. While the transmitter and receiver parameters were substantially the 
same for all three periods there was immediate assurance that the distri- 
bution and frequency meteors would similar the three dates. However, 
evidence accumulated during 1949 and 1950 (12) does suggest that the normal 
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background meteor rates for given hour the day (excluding the effects 
the strong meteor showers) change slowly from day day. would therefore 
reasonable expect that the mean the rates for the periods available 
February and April should representative the rates for March with the 
same antenna system. 


HOURLY RATES 


° 100 140 180 220 260 300 
RANGE KILOMETRES 


Fic. Range distribution curves. 
Solid line—low-gain antenna. 
Dotted line—high-gain antenna. 


The echo rates for the low-gain and high-gain antennae should compared 
only for the range class km. will seen from Fig. dashed curve 
the gain the multiple-element antenna fell off rapidly minimum for 
echoes the range class km. and was appreciably less than the max- 
imum for the range classes km. and 120-140 km. The observed mean 
hourly rates for the km. class were echoes per hour for the high-gain 
antenna and echoes per hour for the low-gain antenna, ratio 1.9 
This figure fair agreement with the measured gain ratio 2.2 for the two 
antennae, view the possible sources error. 


The mean hourly rates, summed over all range classes, for the high-gain and 
low-gain antennae respectively were 307 and 628 echoes per hour. Although 
the high-gain antenna system can detect smaller meteors the low-gain antenna 
covers greater volume the M-region, and hence detects larger total 
meteors despite the decrease echo strength with increasing range. the 
high-gain antenna were tilted 45°, say, the echo rate should increase, both 
because the antenna pattern would intercept larger volume 
meteors coming from radiants high 45° above the horizon could then 
detected, whereas the radiants detectable meteors the overhead zone must 
lie close the horizon. comparing the background echo rates obtained with 
different equipments not only the antenna patterns have specified 
but also the elevation the maximum the antenna beam. Furthermore, 
strong shower present the antenna beam position must completely speci- 
unless the beam lies plane normal the paths the shower meteors 
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very few the shower meteors will detected (6, 7). all further work 
presented this paper the single dipole antenna was used, and the conclusions 
reached will apply specifically such antenna only. 


Variation Rates with Transmitter Power 


series tests was performed ascertain the dependence echo rates 
transmitter peak power. For this purpose the d-c. voltage the power supply 
the modulator was varied discrete steps, and the relative radiofrequency 
power outputs were measured means r.f. current meter the end 
quarter-wave stub line connected across the transmission line the antenna. 
The relative powers could thus measured with accuracy better. 
The absolute powers were known accuracy about 20% only, because 
the computation absolute power required knowledge the pulse width 
and shape, the transmission line surge impedance, and the pulse repetition 
frequency. other parameters the system were changed during the experi- 
ments, except that the receiver was kept tuned the transmitter frequency 
(when the power was varied much the transmitter shifted fre- 
quency 0.1 Mc.). During series observations the transmitter was oper- 
ated high and low power alternatively for half-hour periods and the echo 
counts averaged over several hours. Two experiments were carried out, one 
with power ratio 6.25 (current ratio 2.5 1), see Table IV, 
and the other with power ratio (current ratio 1), see Table 
The echo counts were made for each 20-km. range class from km. 260 km. 
(80 denotes the class 80-100 km., etc.) and the echo ratios are given 
for each range class well for the totals. 


will seen that the ratios total echo counts Tables 
and are reasonable agreement with the transmitter r.f. current ratios, or, 
that the echo rates are proportional the square root the transmitter power. 
will noticed that the values the ratios for the individual km. class 
intervals not show any definite trend from km. 260 km. tenta- 
tively adopt Lovell’s formula (Equation 1b) and assume from the calculated 
antenna pattern that the antenna gain the elevation angle the 
km. range class about one-tenth the gain km. overhead, obtain 


TABLE 
HIGH-LOW POWER TEST 


1900-0415 E.S.T. 30-31, 1950. RATIO 2.5 
HIGH POWER 225,000 Low 36,000 


Range class, km. 100 120 140 160 180 200 220 240 Total 


Hourly rates, high 
power 27. 

Hourly rates, low 
power 

ratio 


% 
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TABLE 


HIGH-LOW POWER TEST 


0000-0400 E.S.T. Nov. 29, 1950. RATIO 4.0 
HIGH POWER 155,000 Low 9700 


100 120 140 160 180 200 220 240 Total 


Range class, km. 


Hourly rates, high 


Hourly rates, low 


theoretical ratio approximately for the minimum detectable electron 
densities the two ranges. This would equivalent range four visual 
magnitudes, accept Herlofson’s proposed relationship. The significance 
the comparatively constant ratios Tables and the range 
increased that the law numbers versus size the meteors the long 
range classes must the same that for the meteors the short range classes 
which are several magnitudes fainter the average. 


Fig. are plotted total hourly echo rates against the square root the 
transmitter power watts. The circled points were obtained from Table 
and the crosses from Table For each the cases straight line may 
drawn through the points and the origin. The lines have not the same slope 
because the data from Table include observations made before midnight 
when the echo rates are lower. the hours after midnight only are used from 
the October period the high and low hourly rates are found 856 and 334 
respectively, and these selected values are plotted crossed circles which lie 
fairly close the line drawn through the November data. This implies that 
the normal echo rate for given hour the day changes but slowly from day 
day, the absence strong meteor shower; conclusion which has been 
borne out several hundred hours observation extending over months, 
made with continuous-wave radio technique (12). error the parameters 
involved computing the absolute power units used Fig. will cause the 
abscissa scale expand contract without changing the relative ratios the 
plotted points. 


The large transmitter could not operated successfully for more than 
variation power. Accordingly, small radar transmitter was built 
determine the ratio echo rates for large difference power, particularly 
compare the echo rate ratio obtained during shower period with that 
obtained nonshower, normal period. During the observing periods the 
low-power radar was maintained 100 output 9.52 m., and the high- 
power radar 9.22 was held closely possible 200 kw. output. 
Separate and identical receiving antennae, receivers, and displays were em- 
ployed the arrangement used for the crossed-polarization experiments, except 
that all the dipoles were oriented N/S. The db. rejection ratio the two 
receivers was thought sufficient ensure that echoes from the high- 
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power transmitter would appear the low-power display, but additional 
precaution the low-power transmitter was triggered 133 after the high- 
power transmitter order delay all low-power echoes km. Under 
normal, nonshower conditions the high-power system should detect 
times many echoes the low-power system (see next 
section for figure-of-merit formula). 


HOURLY RATES 


100 200 300 400 
Fic. Observed hourly rates meteor echoes versus square root transmitter power. 
Circled points—period 1900-0415 E.S.T., Oct. 1950. 


Crosses—period 0000-0400 E.S.T., Nov. 29, 1950. 
Circled crosses—period 0000-0415 E.S.T., Oct. 1950. 


From 7.5 hr. film records obtained between 1700 E.S.T. Jan. and 0200 
E.S.T. Jan. 1951, the observed echo counts were follows: high power 


2287 echoes, low power echoes; ratio This good agreement with 
the predicted ratio for nonshower conditions. 


Jan. 1951, the Quadrantid meteor shower was active. From hr. 
film records secured between 0900 E.S.T. and 2400 E.S.T. Jan. 1951, the 
observed echo counts were follows: high power 6973 echoes, low power 277 
echoes; ratio 25.2 During the period 1200 E.S.T. 1900 E.S.T., when 
the Quadrantid radiant was low the sky, the ratio averaged 19.6 
the radiant rose towards the ‘‘forbidden (12) around the zenith 
fewer shower meteors were detected either radar, and the ratio 
increased. Comparing the shower and nonshower high-low ratios follows 
that the low-power set detected more than twice many echoes during the 
shower would predicted. The hourly rates are themselves 
about twice the normal nonshower rates for the same hours December and 
January, which would make the low-power hourly rates during the shower 
about four times the normal rates. The complete observations the Quad- 
rantid shower, including Doppler velocity measurements, will presented 
elsewhere. 

Variation Rates with Wave Length 


The variation echo rates over wide range wave lengths more difficult 
evaluate than the variation with other parameters. not practical 
attempt rapid change wave length the radar system factor two, 
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say, and the observations are made simultaneously two wave lengths, 
completely separate radar systems are required and the transmitter power, 
antenna coverage, transmission line losses, receiver band width and noise 
factor, and galactic noise should all measured for each system. 
usually difficult make each these parameters identical for the two systems 
(especially the galactic noise factor which varies widely with wave length and 
not under the control the experimenter) some method putting the 
observations comparative basis must used. the same type antenna 
used both wave lengths the evidence the preceding sections suggests 
that the echo rates may equated, far transmitter power and received 
signal power are concerned, the use relative figure-of-merit, 
Here and are the minimum detectable signal powers and are 
assumed include factors such receiver band width, receiver noise, and 
galactic noise. The minimum values quoted Table were determined 


TABLE 
COMPUTED FIGURE-OF-MERIT 
2.0 4.5 5.9 
8.2 


under operational conditions switching the receiver input alternatively 
between the antenna and diode noise generator and equalizing the outputs 
A-scan oscilloscope. The measured values are about db. more above 
the minimum signal that could actually read the film record the 
intensity-modulated display, owing the integrating characteristics the 
film. The measured values will used here because they furnish reliable com- 
parisons between receivers. This expression for consistent with both the 
echoing formulae Equations (1) and (2). The antenna gain factor has been 
left out this discussion because was expedient use similar half-wave 
dipoles all wave lengths, properly scaled size and all oriented 
South. this experiment was necessary measure the absolute power 
each transmitter, hence the error the relative power ratios was the order 


20%. 


The 9.22 radar system was first compared with lower-powered equip- 
ment operating 2.83 The short-wave system was later modified operate 
5.35 m., primarily for some auroral research work. While the equipment 
was available advantage was taken the opportunity compare 9.22 with 
5.35 m., using half-wave dipole each case. The minimum noise powers for 
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the 9.22 and 5.35 receivers were approximately equal, the decrease 
galactic noise (14) the shorter wave length being counterbalanced the 
greater receiver band width. new receiver was built for 5.35 m., but the 2.83 
receiver was early wartime construction and its performance inferior 
similar receivers modern design. The measurements are summarized 
Table VI. addition, the last two columns Table show the figure-of- 
merit values multiplied respectively the three-halves power and the square 
the wave length ratios. 


The observational results are presented Table VII and Table VIII. 
Echoes were counted all ranges the maximum 260 km. available 
the displays. Almost without exception every echo that occurred the shorter 
wave lengths was detected 9.22 The counts and the ratio are given for 
each hour show the consistency the data. 


TABLE VII 
OBSERVED ECHO RATIOS 
Aug. 12, 1950 No. echoes 
Hour E.S.T. Ratio 
9.22 5.35 

0100-0200 1,103 212 5.2 

0200-0300 1,044 234 4.5 

0300-0400 840 184 4.6 

0400-0500 667 131 5.1 

Mean 937 199 4.7 

TABLE VIII 
OBSERVED ECHO RATIOS 
Apr. 3-4, 1950 No. echoes 
Hour E.S.T. Ratio 
2.83 

2230-2330 514 103 

2330-0030 565 

0030-0130 707 101 

0130-0230 691 

0230-0330 654 164 

Mean 626 104 


From Table VII the mean ratio echo counts, 4.7, agrees more closely with 
the predicted ratio 4.5, based Lovell’s formula,* than with the value 
5.9 deduced from the assumption that the minimum echo varies 

manuscript privately communicated Lovell proposed the use sensitivity factor, 
with 100 km. standard, for comparison various equipments, but deleted 


this remark from the published paper (17). The quantity used here expresses the same 
relation relative form. 
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the fourth power the wave length. Comparing the observed ratio Table 
VIII with the theoretical ratios for this case the interpretation indefinite 
both theoretical predictions are below the observed value. reason for the 
discrepancy will suggested the next section. 


Referring Table VII should remarked that the echo counts for the 
hours 0000-0300 E.S.T. include many Perseid meteors normal background 
meteors. 0500 E.S.T. the Perseid radiant was high the sky that only 
the brighter Perseid meteors were being detected. The April period, Table 
VIII, was considered reasonably free any pronounced meteor shower. 


Confirmation Lovell’s Formula 


have examined statistically the possibility that the meteor trail may 
highly efficient the longer wave lengths and have rejected 
this hypothesis and its implied explanation the unusually large numbers 
meteor echoes reported the longer wave lengths. Nor consider that 
multiple echoes from the brighter meteors can begin account for the observed 
rates because the complex meteor echoes comprise only small percentage 
the total number echoes observed the longer wave lengths. Furthermore, 
such multiple echoes are usually clearly recognizable being due one meteor 
and are counted one echo. 


The observed variations echo rates with transmitter power and antenna 
gain suggest direct proportionality between echo rates and the minimum 
detectable value deduced from Lovell’s formula, Equation These 
observations are equally accord with Equation (2), based mechanism 
total reflection from ionized cloud with finite scattering area. However 
the time factor suggests that Lovell’s formula preferred because most 
the small echoes have been shown occur simultaneously with the passage 
the meteor across the line sight (12). The diffraction echo from line 
source detectable without delay, whereas appreciable time would elapse 
before the trail could expand volume having cross-sectional area suffi- 
cient yield echo. Furthermore, the line densities most the detectable 
meteor trails are the order 10° electrons per cm. (assuming 
formula) and likely that the volume density the cloud would drop below 
the critical density before the minimum cross-sectional area was attained. 


The observed variations echo rates with wave length, though not con- 
clusive the other observations, indicate that the rates increase with wave 
length more rapidly than would expected the basis Lovell’s formula 
alone. The explanation may lie the decreasing aspect sensitivity the 
average meteor trail the wave length increased, effect that was strik- 
ingly apparent this work and which has been pointed out Lovell. has 
been shown (12) that the received power intensity the echo diffracted 
from the meteor trail being produced increases proportionally 
the final value given Equation (la), where the distance the 
meteor has still travel along its path before reaching the point (the point 
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where the path normal the line sight). Lovell (8) has put forward ample 
experimental verification that the maximum intensities the echoes from 
given meteor observed two wave lengths will vary \*. However, 
given distance prior the point the ratio the echo intensities will vary 
additional factor other words, the scattering re-radiation pattern 
the meteor trail narrower the shorter wave lengths, hence there 
smaller probability that detectable echo will received. wave lengths 
the range 20-100 m., say, the re-radiation patterns are much broader and 
the effective statistical variation echo power with wave length would 
multiplied that hence the echo rates would vary 


Distribution Nonshower Meteors 


From preliminary analysis 1938 visual nonshower meteors observed 
between 2300 and 0100 E.S.T. the clear moonless nights the Ottawa 
meteor program, Millman has determined the hourly rates meteors var- 
ious magnitude classes. Fig. each solid point represents the total number 
meteors per hour, seen standard team six observers, down and 
including the magnitude given the abscissa the point, the data being 
available half-magnitudes. The visual rates begin fall off the 2nd 
magnitude, for various physiological reasons. solid curve drawn through 
these observational points, terminating the 5th magnitude which the 
limit the visual observers. The solid straight line, drawn through the linear 
part the data, may extrapolated fainter magnitudes represent the 
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Fic. Total hourly meteor rates versus limiting magnitude. 
Solid points—visual observations nonshower meteors. 
Solid straight line—theoretical graph extrapolated from visual nonshower data. 
Circled points and dashed line—graph for 9.22 radar computed from visual-radar 
correlations. 
Dotted curve—proposed theoretical graph for shower meteors. 
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theoretical rates available meteors. The slope the solid line Fig. 2.7. 
Watson (18) mentions the value 2.5 for the increase numbers with each 
magnitude. assumed that there constant ratio numbers any 
two consecutive magnitude classes then the ratio the number meteors 
magnitude the number magnitude the same the ratio the 
sum all meteors from magnitude minus infinity down the sum 
all meteors down The sums have been plotted Fig. for convenience, 
but graph numbers each magnitude versus magnitude will have the same 
slope, though displaced downward factor 1.59 for the slope 2.7 and 1.67 
for the slope 2.5. 


The brightest meteors will detected the radar regardless the path 
orientations. Because the aspect effect the radar will miss increasing 
proportion meteors the magnitudes become fainter. Assume that there 
are meteors per hour given magnitude, theoretically available. this 
number, the visual observers record known number and the radar 
detects unknown number The sampling processes differ, course, for 
the visual and radar meteors but for simplicity let assume that each process 
random. Suppose that visual meteors are definitely correlated 
with radar echoes, methods outlined elsewhere (11). The proportion 
radar meteors the visual observations then assumed the same the 
proportion all radar meteors the total number available, i.e., p/n m/N. 
For meteors brighter than 2nd magnitude, The ratios, p/n, can easily 
determined the 5th magnitude, though the ratios will tend high 
because the increasing probability chance correlation the fainter mag- 
nitude classes are considered. The circled points Fig. show the radar rates 
computed this basis from 607 visual-radar correlations. The dashed straight 
line through the circled points Fig. indicates the tentative relationship 
between rates and magnitudes for the 9.22 radar system operating 
normal power 200 kw. 


Typical rates radar meteors per hour midnight nonshower 
nights would correspond limiting magnitude 7th 8th magnitude for 
the radar system, read directly from the dashed line Fig. For several 
reasons this regarded underestimate the limiting sensitivity 
the radar system. First, the computed points 3rd 5th magnitude 
the radar graph will tend high. Second, the slope the theoretical 
visual line were changed from 2.7 2.5 the corresponding radar line would 
yield limiting magnitude between the 8th and 9th. Third, the radar rates 
were fall off manner analogous the visual, suggested the dashed 
curve Fig. the limiting magnitude would even fainter. The reasons for 
such decrease the radar rates the fainter magnitude would quite 


from the reasons for the visual rate decrease: one suggestion would 


that radar meteors the 9th 10th magnitude ranges tend occur high 
elevations, the upper fringe the M-region (11), where less ionization 
produced for given mass and velocity the particle. Fourth, important 
consideration that has not yet been taken into account the argument far 
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the fact that the visual data include some meteors beyond the range the 
radar gear and many meteors near the limits the radar range where the radar 
sensitivity four five magnitudes less than the zenith. Magnitude cor- 
rections should applied the more distant visual meteors, well, but these 
will smaller than the radar corrections. Expressed other words, the atten- 
uation with range and zenith angle less for the visual observers than for the 
radar. The average range ‘the radar echoes about 130 km. where the sensi- 
tivity about one magnitude less than the zenith. Hence, the radar graph 
Fig. will apply the average radar meteor, and the limit for overhead 
meteors will one magnitude fainter. 


100,000, 


TOTAL HOURLY RATES 


8 


Fic. hourly meteor rates versus limiting magnitude, and versus the nuniber elec- 


trons produced per cm. meteor path. 
Dashed line—graph for 9.22 radar (taken from Fig. 4). 
Solid for 2.83 radar. 
Dotted curve—proposed theoretical graph for shower meteors detected the 9.22 


radar. 


The dashed line Fig. has been transferred Fig. and represents the 
uncorrected relation between normal, nonshower radar echo rates and mag- 
nitudes for the 9.22 system. From Table the theoretical rate ratio for 
9.22 versus 2.83 was computed about whereas the observed 
ratio was about 100 (see Table VIII). suggested above the discrepancy 
may due part the effectively sharper radiation pattern the meteor 
trail the shorter wave length. This would imply that the short wave length 
radar system should lie graph smaller slope than the long wave length 
system. From the 9.22 graph rate 600 per hour corresponds limiting 
sensitivity 7.5 magnitude. the power the 9.22 radar had been lowered 
that the rate was per hour the limiting magnitude would However, 
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the observed rate using the 2.83 radar was per hour, that the graph for 
2.83 should pass through this point the ordinate magnitude 
may assume that the short wave radar will see all the extremely bright 
meteors, and draw the solid straight line, shown Fig. indicative the 
performance radar system 2.83 The short, heavy vertical strokes 
the graphs represent the minimum limiting magnitudes the two equip- 
ments actually used. The 5.35 radar would graph below but close 
the 9.22 line. The limiting magnitude for the 5.35 system estimated 
from Tables and VII about the 6th magnitude and marked 
Fig. with cross. 


The predicted rate based published data the parameters his 
4.2 radar system with half-wave dipole antenna fair agreement with 
his observed nonshower rates, and suggest limit 5th magnitude for his 
equipment. Lovell (1) has stated that his gear has limiting sensitivity 
about 6th magnitude, presumably based the visual scale used Prentice 
and Porter. Clegg and Davidson (3) present data which suggest that the 
English visual observers record large numbers 6th magnitude meteors, 
whereas our visual records contain none. meteors described 5th magnitude 
the Canadian team were considered 6th magnitude meteors the 
English observers then the magnitude scale Figs. and should shifted 
one magnitude the left. Or, slope 2.5 used represent the value 
generally quoted the literature, instead the tentative value 2.7 based 
the observations made Ottawa, similar displacement the limiting sensi- 
tivities towards the fainter magnitudes will obtained, although the dis- 
placement will greater for the fainter magnitudes. sum up, appears 
that the limiting magnitude the 9.22 radar, operating 200 kw. power 
into half-wave dipole, certainly fainter than the 8th magnitude and 
probably the order 9th 10th magnitude. 


The limiting magnitude the continuous-wave Doppler system employed 
Ottawa for the measurement meteor velocities (12) may deduced 
relative basis. The Doppler wave length 9.99 m., the transmitter power 
kw. and the receiver band width ke. per sec. The band width the radar 
receiver 200 kc. per sec. Theoretically, the minimum noise powers the 
two receivers should vary the band width ratio, 100. 
This ratio exactly counterbalances the transmitter power ratio, hence the 
relative figure-of-merit Doppler versus radar should depend only the 
wave length ratio 1.13. The Doppler should thus slightly 
more sensitive than the radar, but the actual Doppler echo counts are found 
vary from 1/10 1/3 the radar echo counts. The discrepancy between 
theory and practice explained part additional source noise being 


-present the Doppler system, owing scattered transmitter radiation 


random phase reaching the receiver reflections from ground objects the 
vicinity. Tree branches waving the wind can increase the Doppler receiver 
noise figure significantly. This effect does not increase the radar noise figure, 
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except for the first few kilometers the range sweep where meteor echoes 
occur. Using the experimental data rather than the theoretical values, fol- 
lows that the Doppler sensitivity less than that the radar one two 
magnitudes, depending local conditions. 


Distribution Shower Meteors 


Evidence was put forward previous paper (11) which suggested that 
the numbers/magnitude relation for strong visual shower does not follow 
the same law for the normal background meteors. The maximum ratio 
shower echoes background echoes during the Geminid shower 1947 was 
also pointed out that the percentage radar echoes duration sec. and 
longer during the Geminid shower was twice great the percentage during 
typical nonshower periods. The conclusion was that there smaller percentage 
faint meteors the shower stream than the normal background, the 
implication being that some agency, such the Poynting—Robertson effect, 
has weeded out the smaller shower meteors over the course centuries. The 
evidence the high-low power experiments with the 100-watt and 200 kw. 
radar system showed that the echo rates obtained during the Quadrantid 
shower with the low-power transmitter were twice the rates that would 
expected the distribution the shower meteors follow the normal 
distribution law. From the predicted nonshower ratio, the observed 
ratio, the limiting magnitude the 100 set estimated from Fig. 
about 3rd 4th magnitude, which well within the visual range. 


Let return and draw the dotted curve tentative represen- 
tation the theoretical numbers/magnitude distribution shower meteors 
well established shower such the Perseids Geminids. the range 
magnitude, the indicated shower rates are three times the normal 
rates, but fainter magnitudes the shower curve arbitrarily drawn ap- 
proach the normal line and cross about the 9th magnitude. Over the 
visual range the ratio shower nonshower meteors would the order 
though there might slight decrease noticed 5th 6th mag- 
nitude. However, 9th magnitude the ratio would about similar 
shower curve for radar meteors drawn the dotted line Fig. again 
crossing over about the 9th magnitude. The radar rates the 200 kw. 
9.22 equipment would then have doubled during strong shower while the 
visual rates would have quadrupled, qualitative agreement with observation. 
Radar equipment having lower sensitivity, whether virtue lower power 
shorter wave length, would record increase echo rates during showers 
more proportion the increase visual meteor rates. Lovell has stated 
that his equipment 4.2 does indeed yield echo rates proportion the 
visual rates. 


Olivier (15) has concluded from his studies telescopic meteors that the 
meteor rates begin depart from the expected law increase number 
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the 7th 8th magnitude, and fall off rapidly after the 9th 10th magnitude. 
does not state explicitly what fraction his observations might shower 
meteors: the proportion were large his findings would qualitative agree- 
ment with our own. However, most his observations are nonshower 
meteors would imply that our theoretical numbers/magnitude relation for 
normal nonshower meteors (solid line, Fig. 4), and also our normal radar graph 
(dashed line, Fig. Fig. 5), should commence curve downwards the 
7th 8th magnitude. this were the limiting sensitivity the high-power 
radar would even fainter than estimated above, possibly the order 12th 
magnitude, but the high-low power observations not lend support the 
suggestion. 


interest speculate upon the performance kw. radar set 
operating 100 m., say. Despite the lower power the echo rates should 
increased factor over the 9.22 gear. The 100 set would oper- 
ating line above that for the 9.22 set and below the theoretical graph 
available meteors. Roughly speaking, the sensitivity the low-power 100 
set would more than two magnitudes fainter than the high-power 9.22 
set, and the over-all meteor rates seen during well known showers the long 
wave set would show little increase over normal rates, although the per- 
centage long-enduring echoes would increase. 


the other hand, the earth were pass through clouds meteoric dust 
one might expect very great numbers extremely faint meteors occur, which 
would escape detection high-power radars shorter wave lengths 
but would seen long wave equipments. The echo rates would high 
that the resolution the conventional ionospheric display would inadequate 
show individual echoes, and continuous band echo should obtained. 
This purely speculative but may offer explanation for some sporadic-E 
type echoes detected ionosphere equipments moderate power working 
the range 50-150 any sporadic-E type echoes were due these dust 
clouds the echoes should appear simultaneously over hemisphere the earth. 
Pineo (16) has indicated that sporadic-E echoes general cannot correlated 
with the well established visual showers. Our own observational data confirm 
his conclusion, which compatible with the argument advanced above that 
the visual showers should have little effect ionosphere equipment. The sug- 
gestion that showers faint meteors can occur which would affect ionospheric 
records may worth some investigation. These would not necessarily 
meteor showers having predictable orbits the usual sense. They would more 
likely due chance encounter with relatively dense clouds dust, possibly 
interstellar material. 


Using Equation found that the minimum line density electrons 
that can detected the 9.22 radar the order 10° elec- 
trons per cm. may associate the value 10° with the 9th magnitude, 
estimated above the limit for the 9.22 set. was shown that meteor 
rates varied the square root the power, directly with hence the 
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complete correlation between and visual magnitude may easily estab- 
lished. This relation shown the second abscissa scale Fig. inde- 
pendent line reasoning led Herlofson (5) estimate that 10” would 
correspond Ist magnitude meteor and toa 6th magnitude meteor, 
which agreement with the scale Fig. 


Summing up, has been demonstrated that Lovell’s formula accord 
with the experimental evidence over wave length range 2.8 m., when 
consideration given the effective broadening the radiation pattern 
the meteor trail the wave length increased. difficulty encountered 
explaining the higher meteor rates observed the longer wave lengths, 
and there appears need postulate additional echoing mechanisms. 
particular, the mechanism based Herlofson’s plasma 
theory has received support. The increase visual rates over radar rates 
during strong visual shower, compared the rates observed during normal 
nonshower periods, can explained assuming that there are fewer small 
meteors the shower stream. The argument has tentatively been extended 
wave lengths 100 suggest that echoes might obtained from 
clouds meteoric dust which could detected only long wave radars, 
with the appearance sporadic-E type echo. 
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Abstract 
shown that the equation 


where pressure function and resistive frictional force, taken 
the fundamental dynamical equation, possible explain diffusion solu- 
tions, solid solutions, and adsorbing solid well the interdiffusion 
diffusion both liquid and solid solutions the pressure function 
the osmotic pressure; diffusion adsorbed gases the pressure function 
the spreading pressure the adsorbed layer. Since these pressure functions 
are thermodynamic properties the system the pressure gradient independent 
the physical picture assumed for the system. evaluate Ax, however, 
necessary assume specific physical mechanism for the movement the 


molecules. 


Recently (1) have shown that if, place law, giving the flux 
diffusion proportional the concentration gradient, assumed that 
pressure gradient the operative force for diffusion, then possible obtain 
reasonable picture the diffusion gases through solids. that time 
were mainly interested finding some explanation the flow adsorbed 
gases through adsorbent solids. Since that paper was written, have become 
through recent publication (4), the work that has been done the 
diffusion metals metals and the theories that have been put forward 
explain such diffusion. After consideration this work, became apparent 
that extension the ideas presented earlier would lead the same equa- 
tions those adopted explain the diffusion metals and would, conse- 
quence, present single picture covering all aspects diffusion. shall, 
then, this paper, apply the fundamental equation developed earlier dif- 
fusion liquid and solid solutions and show that this equation leads the 
same expression the current theories such diffusion. 


Fundamental Equation and Application Gases 


the earlier work was shown that derive the differential equation for 
diffusion was necessary combine 
(i) dynamical equation defining the motion, 
(ii) equation continuity expressing the law conservation 
mass, 
and (iii) equation state. 


the past Fick’s law 
(1) 


received September 18, 1950. 
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has been implicitly used the fundamental dynamical equation. Here 
the rate permeation the flux diffusing, the diffusion coefficient and 
the concentration. From hydrodynamical considerations the fundamental 
equation should have the form 


z= 0, (2) 


where potential function having the properties pressure and 
the resistive frictional force that opposes the motion. This equation has 
general validity for diffusion when inertial effects and external forces can 
neglected whereas law secondary law derived from under the 
special conditions arising when directly proportional apply Equa- 
tion (2) the various types diffusion necessary only choose correctly 
the pressure function that operative and express the resistive forces 
terms the physical properties the system. 


was shown long ago Maxwell (5) that the application Equation (2) 
the interdiffusion two gases presents difficulty and, fact, leads 
directly Fick’s law when the gases obey the gas law. Maxwell’s fundamental 
equations for the interdiffusion two gases are 

where pip2, and are the partial pressures, the densities, and the 
velocities the two gases and respectively and the coefficient 


resistance. represents the force that must applied unit mass the gas 
maintain unit velocity when the density the second gas unity. 


easy prove that 


Cp2P Ox Cp, P Ox ( 
assume that the gases obey the gas law that 
2 2 
place Equations (6) are equivalent Fick’s law. Thus, 


the gases behave perfect gases, Fick’s law valid expression for the 
diffusion but when the laws perfect gases are not applicable necessary 
Equations (5) for the fundamental equation. 


wae 
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Diffusion Solutions 


may apply the same concepts that were used above for the interdiffusion 
start with the fundamental equation 


the theory solutions has long been known that the osmotic pressure 
measure the tendency solute dilute itself and early 1888 
Nernst (6) assumed that diffusion solution was proportional the gradient 
the osmotic pressure. similar approach has been used Sackur (8), and 
Bryan and Porter their article diffusion the 14th edition the En- 
cyclopaedia Britannica also suggest that solutions the osmotic pressure 
the operative force for diffusion. is, therefore, not only reasonable 
assumption but one that has considerable support assume that the osmotic 
pressure the pressure function that should used for solutions Equation 
(2). Before replacing the partial pressures gaseous diffusion with the osmotic 
pressures solutions necessaty, however, consider carefully the way 
which the term osmotic pressure used. 


The osmotic pressure has been defined Guggenheim (2) follows. 
solvent and solution are separated membrane that allows the solvent 
pass without hindrance, but prevents the passage the dissolved sub- 
stance, the solvent will tend flow through the membrane into the solution. 
The hydrostatic pressure that must applied the solution stop the flow 
solvent through the membrane the osmotic pressure the solvent. 
the difference the pressure the solvent the pure phase and the 
solution. obvious that the pressure tending make the solvent mol- 
ecules diffuse equal the osmotic pressure. also apparent that, 
could find semi-permeable membrane that would allow the passage the 
solute molecules but prevent the passage the solvent molecules, exactly the 
same considerations would hold and the osmotic pressure the solute atoms 
would measure the tendency the solute molecules 
this case the pressure responsible for the diffusion would equal the 
pressure solute molecules solution since their pressure the pure solvent 
would nil. 


this picture, the osmotic pressure defined the hydrostatic pressure 
that must applied stop diffusion; follows, therefore, that the pressure 
causing the diffusion would equal but opposite sign the osmotic pres- 
sure. interesting note that, since the concept semipermeable mem- 


brane could applied equally well mixture gases, the partial pressures 


would play there the role assigned solutions the osmotic pressures. 


determine proceed for gases and assume that are the mass 
solute and solvent respectively per cc., and their velocities and.C the 
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coefficient resistance, then the portion the solvent contained small 
element volume will experience from the solute resistance 
and the solute contained the same volume element will experience from the 
solvent equal and opposite resistance With these consider- 
ations write for the equations flow 


where the osmotic pressure defined above. The positive sign used since 


the resistive force the same direction the increase pressure. 


solution the hydrostatic pressure will remain constant that the net 
transfer mass across plane must zero. 


where the total mass per cc., that is, the density the solution. 
Substituting (7), get 
Transposing obtain our equation diffusion giving the mass flux diffusing 
across plane 


x 


holds then can write 


and our equations become 
RT Op. Ope 


Thus see that, for solutions that obey van’t Hoff’s law, Fick’s law valid. 
noted, however, that Equations (11) are generally applicable, 
whereas Equations (13) are restricted ideal solutions. find general 
relation between diffusion and concentration gradient necessary intro- 
duce the relation between and the concentration for nonideal solutions and 
that shall now proceed do. 


Equation Diffusion Terms Concentration Gradient 


According Guggenheim (2) the osmotic pressure the solvent non- 
ideal solution given 
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where the partial molar volume the solvent, the partial molar 

fraction, and the activity coefficient. must remembered that the 

activating pressure for diffusion opposite and equal this. Therefore 

can write 


Here the partial molar volume solute. 
The chemical potentials the two substances are given 
where and are constants depending only temperature and pressure. 
From (15) and (16) have 


Substituting Equation (11) get 


= and pot, = — 2H? (19) 


see, therefore, that the mass flux diffusion proportional the 
gradient the chemical potential. 


Several workers have suggested that diffusion should taken propor- 
tional the gradient chemical potential rather than the concentration 
gradient (cf. Onsager and Fuoss (7), Jost (3), LeClaire (4)). Although this move 
seems have been made intuitively these workers, the above analysis 
shows that reasonable deduction from first principles. Proceeding, 
have from (16) 


number solvent molecules and total number molecules per 
cc. Hence 
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may consider constant that 


Similarly 


Substituting (18) get 


These equations are equivalent those given Onsager and Fuoss (7). 


Introduction the Mobilities 


Since the past has been customary use the concept mobility 
treating diffusion solutions, interesting introduce the mobility into 
our equations. define the mobility the velocity flow under unit force. 
Since our fundamental assumption that the gradient osmotic pressure 


the force the particles, the mobility obtained putting Equa- 
tions (11). 


(25) 


Substituting (23) obtain our equations diffusion terms the 
mobilities 


and 


obtaining these formulas have defined the mobilities according the 
picture diffusion that has been developed here. This, course, not the 
procedure that has been used hitherto and, for comparison, interesting 
investigate the form which the equations would take had assumed, has 
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been done the past, that the driving force for the diffusion was the gradient 
chemical potential rather than the osmotic pressure gradient. 


This value for similar that used Onsager and Fuoss (7). 


Proceeding the assumption that defined (27) should get 


= + hi) dor and = 1 fs) (28) 


which the equation given Onsager and Fuoss (7) terms the mobility. 


Equations (26) and (28) differ the factor complete the picture, 


necessary see where this factor arises. Clearly the result the two 
ways which has been defined—in one picture the force has been taken 


pressure gradient represents force tending towards expansion but the 
meaning spatial gradient chemical potential not immediately evident. 
The osmotic pressure measure the force tending make one consti- 
tuent solution change its volume whereas the chemical potential 
measure the force tending towards change concentration. change 
volume will result change concentration and Guggenheim (2) has 
pointed out that bears the same relation that bears Therefore 


reasonable assume that equivalent Since, except for con- 


stant factor, equal one dimensional flow, and since quite 


obviously the force behind any diffusion, follows that rather than 
the true representation terms the chemical potential the force causing 
diffusion. Thus the two equivalent equations are 


and 


That these are equivalent immediately evident from the fundamental ther- 


The origin the now obvious. have 
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since definition the partial molar volume. Since again one dimen- 


sional flow the change independent variable from involves only 
constant factor our result obtained. 


Diffusion Solute into Solvent Assuming that the Solvent 
Stationary 


deriving these equations for diffusion solutions have assumed that 
both solvent and solute diffuse or, what equivalent, that there dis- 
tinction between solvent and solute. dilute solutions negligible error 
introduced assuming that the solvent present much greater quantity 
than the solute and that may considered stationary that the dif- 
fusion confined the solute particles. Since the resulting equations are 
simpler than those produced above and have considerable application 
worthwhile derive them. 

before assume that the driving force for the diffusion but 

now the osmotic pressure the solute. The resistive force may written 
Cpu where now includes term for the density the solvent; the con- 
centration solute per cc. and the velocity solute with reference the 


solvent. Thus 
(29) 


d 30 
(30) 


comes 


before and the equation be- 


(31) 


terms the chemical potentials the same considerations above hold 


that 


and the equation terms the concentration gradient becomes for nonideal 
solutions 
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The mobility now given 


that the reciprocal the resistance per cc. 


The diffusion equation for nonideal solutions terms the mobility 


“= 


Diffusion Metals Metals 


obvious that have solid solution one metal another, can 
apply, with very few changes, all the considerations given above for liquid 
solutions. This is, effect, what has been done the recent treatment 
diffusion metals which has been ably reviewed Claire (4). 
true that solid solutions the concept osmotic pressure has not been de- 
veloped but there can doubt that the molecules the solute metal are 
free change their positions the solvent metal they will exert pressure 
similar osmotic pressure and this pressure will direct measure the 
tendency for diffusion. Such pressure can evaluated thermodynamically. 
metals, however, the chemical potential quantity that more readily 
available, and, consequence, probably preferable express the dif- 
fusion function the gradient chemical potential. Since the gradients 


chemical potential and osmotic pressure differ only the im- 


material which used. The concept osmotic pressure useful, however, 
relating the diffusion coefficient the general hydrodynamic equation 
flow where the resistive forces are balanced pressure gradient. 


Once fundamental dynamical equation for the diffusion metals obtained, 
the point great interest relate the coefficient resistance or, what 
equivalent, the mobility the physical properties the metal and speci- 
fically the different models atomic migration. this paper, however, 
are not concerned with that problem and shall proceed further with it. 


Diffusion Gases through Adsorbing Solids 


earlier paper (1) have shown that Equation (2) may applied 
the diffusion adsorbed gases through solids the spreading pressure the 
adsorbed layer taken the operative force for the diffusion. This idea fol- 
lows logically from the unified picture presented here since the spreading 
pressure adsorption thermodynamically equivalent the osmotic pres- 
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sure solutions. Since the pressure may obtained from the thermodynamic 


relation 
(36) 
OV/T.n 
where the free energy, matters not, determining the pressure gradient, 
whether the relation between the gas and the solid assumed that 
adsorption solution. realize this eliminate immediately much 
the confusion that has existed the past all theories diffusion. 
evaluating however, necessary assume some physical picture for the 
diffusion since only through such picture possible relate the resistive 
force the physical properties the system. 
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THE DIFFUSION ADSORBED GASES 
THROUGH 


Abstract 


The equations for the diffusion adsorbed gas through solid are developed 
the assumption that the resistive force directly proportional the number 
molecules adsorbed, that is, that the resistive force per molecule constant. 
The equations are compared with experiment and also with the equations de- 
veloped earlier paper where the assumption was made that the resistive 
force per molecule inversely proportional the number molecules adsorbed. 
The equations developed here are not satisfactory those developed earlier 
and pointed out that this consequence the fact that heats adsorption 
have been found experimentally decrease the amount adsorbed gas 
increases. 


earlier paper (1) have shown how the equation 


where pressure and resistive frictional force, may applied 
the diffusion adsorbed gases through solids and for such diffusion con- 
cluded that the operative pressure the spreading pressure, the adsorbed 
layer. the second paper (2) the use Equation (1) was extended all 
types diffusion and was shown that, general, the operative pressure 
can determined from the thermodynamic relation 


(2) 
where the free energy and the volume. the use (2) the pressure can 
determined independently any physical picture the relation between 
the diffusing gas and the solid. the other hand, determine the resistive 
force necessary assume some mechanism flow for the diffusing 
substance. If, therefore, Equation (1) accepted the fundamental dynam- 
ical equation for diffusion, the application any specific system resolves itself 
into evaluating correctly terms the known physical properties the 
system. 


the treatment the diffusion adsorbed gas through solid 
originally given (1), was assumed that the resistive force could put equal 
where the coefficient resistance and the resultant velocity 
with which the adsorbed layer moves over the surface the solid. this 
assumption treat which the resistive force per unit volume when the 
velocity unity, constant independent the amount gas adsorbed. 
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was originally pointed out that taking the resistive force independent 
the amount gas adsorbed was assumption that would require experi- 
mental confirmation. The experimental results quoted that paper did indi- 
cate that such assumption was justified and the resulting equations were 
put forward solution the problem. 


expanding the scope Equation (1) cover diffusion liquid and solid 
solutions, has become apparent that, putting —Cu, were not 
making the assumption that was intended. Our intention had been make 
the resistance per molecule, rather than per cc., independent the amount 
gas adsorbed but, owing wrong interpretation the fundamental equation, 
have not obtained that result. this paper, therefore, shall develop 
the equations originally intended and compare them with experiment and 
with the earlier equations. 


the assumption that the resistive force per molecule constant inde- 
pendent the number molecules adsorbed write where 
the coefficient resistance per molecule and the total number 
molecules adsorbed per cc. The fundamental equation becomes 


Hence the number molecules diffusing across unit area unit time given 


shall now apply this equation different models adsorption. 


Diffusion Mobile Monolayer 


Boltzmann’s constant, and the area surface per cc. Substituting (4) 
get 
follows from this that, for mobile monolayer, the diffusion directly pro- 
portional the concentration gradient. Thus the diffusion obeys Fick’s law. 


Since mobile monolayer the two dimensional analogue perfect gas 
not surprising find that Fick’s law applicable. 


(5) 


Diffusion Ideal Localized Monolayer (Langmuir Adsorption) 


For adsorption that can represented Langmuir isotherm the spread- 
ing pressure given the relation 


(6) 


log 


o= 


. 
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where the number sites per cc. and Substituting (4) 
get 

(7) 


(8) 


Comparing with Fick’s law the form where the concen- 


tration and the rate permeation, see that 


(9) 
This the relation between the diffusion coefficient Fick’s law and the 
Langmuir’s equation. 


interesting find the expression for the differential coefficient terms 
the gas pressure. The Langmuir equation gives 


aN N,; Op 


and the expression for the diffusion becomes 


CnA 


Comparing this with diffusion equation the form see 
that 
(11) 


apparent, therefore, that becomes less increases and when 
large compared unity finally approaches zero. Such behavior 
expected from the fact that Langmuir adsorption approaches saturation 
increased. 


practice always integrated value the diffusion that measured. 
order compare theory with experiment, assume that the diffusion 
takes place across barrier thickness when there concentration 
steady flow, the equation continuity, constant. Then 
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N=No 
(12) 


most experiments the pressure one side the barrier zero take 
Therefore 


Nu = log —y 


log apo). 


test this equation compare with the measurements Smithells and 
Ransley (7) for the diffusion hydrogen through copper 450°C. These 
measurements were also used test the equation developed the earlier 


Oo 20 40 60 60 100 120 440 


Fic. Comparison theoretical curve with the measurements Smithells and Ransley for 
the diffusion hydrogen through copper 450°C. Measurements. 


paper. The agreement between theory and experiment shown Fig. 
calculate the curve, has been taken equal 0.025 and the curve has been 
made agree with the experimental results the highest pressure measured, 
that is, 129.5. With these adjustments, the theoretical curves can 
said explain satisfactorily the experimental results but should noted 
that this equation does not give the reversed curvature suggested the 
measurements very low pressures and reproduced the equation developed 
the first paper (1). the basis these measurements seems that the 
original equation slightly more satisfactory. 


pointed out the earlier paper the real test these equations will come 
only when possible measure independently adsorption measure- 


40 ° 
6 
Nu 
4 
. 
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ments; present adsorption isotherms are available the temperatures 
which diffusion measurements are possible. not surprising that both 
equations are satisfactory explanations these experimental results since they 
differ only low pressures where the measurements are hardly definite enough 
discriminate between them. 


Diffusion Gas Obeying the BET Isotherm 


The spreading pressure gas adsorbed according the equation 
Brunauer, Emmett, and Teller (the BET equation) given 


(14) 


where the total number localized sites and the number these 
sites that are occupied. Substituting this value for Equation (4), get 


BkT 
(15) 
CnA 
eliminate use the relation 
where the total number molecules adsorbed and constant. get 
(17) 


which can written 
(18) 


practice the gas pressure invariably measured rather than the con- 
centration and express the diffusion terms the pressure have the 
BET equation the form 


(19) 


s 
(20) 


that the equation diffusion terms the pressure becomes 


q 
— 
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2 
(21) 


Thus write the equation diffusion the form 


have 
RELATIVE HUMIDITY 
Fic. differential coefficient with pressure when the adsorption follows the 


BET isotherm. 


4 
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cBkT 


Fig. have plotted against The curve shows that 


m 


decreases low humidities from its value zero pressure, passes through 
minimum and then rises very sharply the pressure approaches saturation. 
see that low pressures the diffusion similar that obtained for 
Langmuir adsorption and this expected since the BET isotherm low 
pressures follows the Langmuir isotherm. high pressures, however, the 
diffusion increases infinity following the BET isotherm which gives infinite 
adsorption saturation pressures. 


Again, obtain comparison with measurements, necessary use 
integrated form the equation. Putting 


can write Equation (17) the form 


(24) 
BRT 


consider again the diffusion through barrier thickness with 
that steady flow constant. Then 


test this equation make use the measurements King (5) the 
diffusion water vapor through keratin. King gives curve showing the varia- 
tion rate flow water vapor through horn keratin with change relative 
humidity the high pressure side. addition gives the adsorption iso- 
therm for the keratin that possible evaluate and directly from 
the BET equation. 


Fig. the theoretical isotherm the BET equation using 5.66 and 


17.67 compared with the experimental isotherm King. Above 


relative humidity 0.5, the BET theory gives higher adsorption than 
found experimentally but below 0.5 the agreement satisfactory. 


The experimental results King for the diffusion water vapor through 
keratin are reproduced Fig. along with the diffusion calculated according 
Equation (26) (full line curve). The theoretical curve has been made agree 
with the measurements relative humidity 0.5. Although below 0.5 
relative humidity the same general trend can seen theory and observation, 
above this value the calculated results rise much too slowly. Although 
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this region that the BET equation deviates from the experimental isotherm, 
does not seem possible attribute the difference between the theoretical 
curve and the measurements wholly the inadequacy the adsorption theory. 
would appear that the BET theory would make the diffusion higher rather 
than lower the high relative humidities. can only conclude that Equation 
(26) not satisfactory explanation the diffusion water vapor through 
keratin. 

° 


MOISTURE CONTENT (PER CENT) 


RELATIVE HUMIDITY 


Fic. the BET isotherm with the measurements King for the adsorption 
water vapor keratin. Measurements. Curve calculated with 5.66 and 17.67. 


Since the agreement between Equation (26) and the experiments King 
cannot considered satisfactory, interesting revert the equation 
developed the earlier paper. That equation was based the assumption 
that Cu. The equivalent equation (26) has the form 


When the data King for the adsorption water keratin are applied 
this equation obtain the dotted curve Fig. This curve better 
agreement with the experimental results but still gives values for the diffusion 


° 
° 
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that are too low high humidities. However, the fit with measurements 
relative humidities below 0.5 good and the general shape the curve 
follows closely that the measurements that reasonable believe that 
this theory very close the true picture. 


Since for both the Langmuir and the BET adsorption, the earlier theory, 
putting A,= Cu, better agreement with experiment than the theory 


RATE FLOW (KING) 


RELATIVE HUMIDITY 


Fic. Comparison theory with the measurements King for the rate flow water 
vapor through keratin. Measurements. Full curve calculated from equation 


dotted curve calculated from equation 


this implies. Since the resistance per unit volume, the relation 
equivalent the assumption that molecule contributes resis- 
tance the flow that independent the number molecules adsorbed 
whether large small. This the same thing assuming that mole- 
cules adsorbed low pressures are bound the adsorbing surface with 


/ 
/ 
2 4 
7 
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the same energy those molecules adsorbed near saturation. That this not 
fundamental assumption the BET theory and not surprising, 
therefore, that the equation deduced from this assumption has failed agree 
with experiment. 


the other hand, the assumption that equivalent saying 
that the resistance per cc. independent the number adsorbed molecules. 
The contribution made each molecule this picture equal C/N and 
therefore inversely proportional the number molecules adsorbed. This 
amounts the assumption that the binding force between adsorbed molecule 


and absorbing surface varies There have been several estimates the 


variation heat adsorption with number molecules adsorbed for 
instance, Babbitt (3), Guthrie (4), Miller (6)) and all have found that the heat 
adsorption decreases from high value zero adsorption value 
saturation approaching that the heat liquefaction the gas. general, 
the curves relating heat adsorption number molecules adsorbed are 
such shape that inverse relation between heat adsorption and number 
molecules adsorbed would fair approximation these curves. This, 
then, is, without doubt, the reason why such good agreement found for 
diffusion when the assumption made that Cu. the basis Fig. 
however, apparent that BET adsorption the molecules high moisture 
contents are even more mobile than required the relation 1/N. 
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ESTIMATE THE RELATIVE HELIUM CONTENT 
THE EARLY-TYPE 


ANNE UNDERHILL 


Abstract 


The relative abundance hydrogen helium O-type atmospheres 
estimated comparison the computed strength the line with 
atmosphere computed when the ratio 5.67 and when 13.33, and 
the equivalent width the line derived. Observations the strength \4541 
the spectra the star, Lacertae, and the are 
given, and the observed equivalent width \4541 real stellar atmosphere 
spectral type corresponding the O9.5 model atmosphere estimated. 
comparison the observed strength \4541 with the computations 
indicates that the abundance ratio lies the range 25, that is, the 
relative abundance helium reduced factor 1/3 1/2 from the 
amounts estimated several authors other spectroscopic methods. 
shown that the 09.5 model atmosphere gives fair representation late 
O-type main-sequence atmosphere, and that the computations can best 
reconciled with observation reducing the relative abundance helium 
about number, the rest the stellar atmosphere being effectively 
pure hydrogen. The computations for the model atmosphere are not 
agreement with this result, but reasons are given for believing that the model 
atmosphere does not give good representation early O-type stellar atmos- 
phere, and that consequently the results obtained from this model atmosphere 
are not equal weight those obtained from the O9.5 model atmosphere. 


Introduction 


problem general interest derive the relative abundance helium 
hydrogen stellar atmospheres. One way doing this means model 
atmospheres. Once model atmosphere, calculated under certain well defined 
assumptions, has been obtained, the strength any spectral line produced 
may computed function the abundance the carrier atom the 
atmosphere. plot the computed equivalent width the line against 
abundance defines theoretical curve growth for the line. Comparison with 
the observed strength the line real stellar atmosphere corresponding 
spectral type the model gives estimate the abundance the element. 


stellar atmospheres the line \4541 Hell will broadened Stark 
effect the same way are the hydrogen lines. statistical theory the 
Stark broadening the hydrogen lines has been worked out Verweij (16), 
and shall use this theory. Following (13), shall consider only 
the Stark broadening caused protons. Neglect the Stark effect due 
electrons means that the broadening Stark effect the line will 
underestimated and that result the relative abundance helium derived 
will, anything, little too large. 


have reviewed the statistical theory the Stark effect hydrogen 
previous publication (10), and find (10, Equation (15)) for line which 
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the Stark broadening each component greater than the broadening due 
the thermal motion the atoms and radiation damping, that the line- 
absorption coefficient per gram star material corrected for stimulated 


emission, 
2 2 


mc? Avp 
Here the number atoms per gram star material excited the lower 
level the line, the normal field strength e.s.u., the part the 
total f-value the line due the Stark components which are displaced 
electric field, the part the total f-value due components zero- 
displacement, S(a) the function introduced Verweij (16) describe the 
Stark broadening hydrogen being the displacement expressed 
the function (a, gives the shape the line under the action 
thermal broadening and damping dispersion type, and the half- 
width the Doppler core the line. use Equation (1) for \4541 Hell, 
S(a), and must obtained. The function (a, may found from 
the tabulation Harris (6). 


The definition the function S(a) 


(2) 


where the summation carried out over all the Stark components positive 
and negative displacement. Here the f-value the Stark component 


(3) 


where the square the matrix element for the transition question, 
and the Holtzmark distribution function, which tabulated, for 
instance, Verweij (16). The quantity numerical factor depending 
the quantum numbers the Stark component, 


where the nuclear charge, this case, and 


Equation (4), the chief number the upper level the line 
and and are the parabolic quantum numbers particular upper sub- 
level; are similar quantities for the lower level the line. For \4541 


(5) 
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where the summation extended over Stark components positive and 
negative displacement. The quantity defined similar manner, the 
sum being extended only over components zero displacement. 


have computed the necessary matrix elements and evaluated the 
according Equation (3), using the formulae for the matrix elements given 
Gordon (5) and incorporating the correction the misprint previously 
pointed out (10). From this material the function S(a) for \4541 Hell has 
been computed according Equation (2). This function, which given 


TABLE 
THE FUNCTION S(a) For 


0.000 0.14 1.332 108 
0.005 8.105 0.18 1.034 
0.01 1.322 108 7.854 
0.03 1.604 108 0.30 4.174 
0.04 1.593 108 2.420 
0.08 1.754 108 0.80 2.706 
0.10 1.648 108 1.00 1.477 
0.12 1.502 108 1.50 5.00 105 
0.14 1.332 108 2.00 2.11 105 


Table has generally the same shape the similar functions computed 
Verweij (16) for Ha, Hy, and find that 0.018109 and 
0.000539. 


the wave-length range for which the collisional theory line broadening 
valid, that is, within the distance the line center, the line absorption 
coefficient 

where the total f-value the line, and tne value the damping constant 
(which enters into the quantity includes the effect collisions. For the 
late O-type stars, the distance AX, is, following the arguments (14), 
found almost 0.2 Hence may use the statistical theory line 
broadening over effectively the whole line profile. 


The Model Atmosphere Computations 


The method computing line profiles has been described elsewhere 12). 
Using the form the line-absorption coefficient given above, Equations (1) 
and (6), can obtain each level the model atmosphere where and 
are known. The ionization and excitation the helium ions assumed 
take place accordance with state local thermodynamic equilibrium. 
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There are two O-type model atmospheres for which shall compute the 
strength the line \4541. The first (11) has 25,200°K. and log 4.200. 
The integrated flux from this model corresponds effective temperature 
36,800° and the computed relative strengths lines and 
SilV place the spectral type this model The second model atmos- 
phere (12) has 35,000° and log 4.200. The integrated flux from this 
model gives effective temperature 44,600° and the spectral type, according 
the absolute strength the line Hy, about O5. The structure both 
models has been computed using the equation mechanical equilibrium 
(allowance being made for the radiation pressure and both models 
the stellar atmosphere assumed composed hydrogen and helium 
only, the ratio 5.67 number. 


The structure the model atmosphere, that the temperature, electron 
pressure, and total pressure any optical depth depends upon the com- 
position the model atmosphere; hence, strictly, should compute new 
model atmosphere for each value the ratio assumed. However, since 
preliminary calculations showed that the structure was not very sensitive 
the abundance helium assumed, have assumed that the structure the 
model atmospheres (11, 12) found for the case 5.67 gives ade- 
quate representation the structure when the ratio varies somewhat 
from 5.67. 


Table the computed profile the model atmosphere 
given for the case 5.67 and for 13.33. The equivalent 


TABLE 
THE LINE PROFILE THE MODEL ATMOSPHERES 


| 


0.0 0.778 0.0 0.848 844 
0.4 0.872 0.906 0.5 863 0.873 
0.891 0.917 1.0 0.895 
2.0 0.940 0.951 5.0 0.991 0.994 
5.0 0.977 0.981 10.0 0.997 0.997 
10.0 0.992 0.992 


width the first case 0.81 equiv. and the second case 0.65 equiv. 
Also given the computed profile for the model atmosphere for 
these two cases. For the model atmosphere the equivalent width 0.71 
order compute the line-absorption coefficient 0.0 and compute 
the part the line-absorption coefficient due the Stark components zero 
displacement, was assumed, accord with previous results, that the damping 
constant, was times the classical damping constant. Since the exact value 
used for the damping constant immaterial, was not considered worthwhile 


| 
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evaluate this constant detail; instead value the expected correct order 


The Relative Abundance Helium Hydrogen 


The computed equivalent widths for the model atmosphere are 
plotted Fig. give curve growth. The ratio 5.67, found 
(15) for the star Scorpii and also the value 10, found 
Aller (1) for Lacertae and Bowen and Wyse (4) and Aller and Menzel 
(2) from study the relative intensities the hydrogen and helium emission 
lines planetary nebulae, are indicated. The curve growth has been drawn 
straight line through the two computed points. 


MODEL ATMOSPHERE WITH 


25,200° 
LOG 


units equivalent width, are equivalent angstroms. The relative abundance 


number. 


shall now obtain estimate the observed equivalent width \4541 
spectra and, interpolation Fig. derive the abundance 
helium late O-type stellar atmospheres. 


The equivalent width \4541 has been measured five plates 
Lacertae, the measurements being made integrating the area 
the line intensity tracings the spectral region surrounding \4541. The 
results these measurements and measurements four plates Orionis, 
are given Table The prism spectrograph used and the dispersion 
\4541 each plate are indicated. Each measure weighted according 
the dispersion and the quality the spectrogram. The weighted mean equiv- 
alent width given the bottom Table III for each star. 


For Lacertae have found the following published measures the 
equivalent width \4541. Weights have been assigned these measures 


w 
0-9 
0-8 
0-7 
0-6 
0-5 
0-4 
0-0 
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Author Weight 
Aller (1) 0.37 equiv. 
Petrie (9) 0.68 
Baker (3) 0.79 
Weighted mean 0.63 


according their seeming reliability and weighted mean equivalent width 
has been obtained. This value, 0.63 equiv. rather smaller than the value 
derived from the new measures Table III, and probably this difference 
may accounted for, least part, the neglect the extended, shallow 
Stark wings \4541 (see below) the above authors, although the method 
adopted Baker should display these wings. The discrepancy upsetting, 
but study our high-dispersion line profiles, and the computed shape 
the line \4541, encourages give considerable weight our measure- 
ments. The weighted mean equivalent width Lacertae from all 
the measures 0.85 equiv. The line \4541 not present (15) the 
star, Scorpii, hence likely that the equivalent width \4541 star 


TABLE III 
THE MEASURED EQUIVALENT WIDTH 


Lacertae, O9V Orionis, O9.5V 
| 

Pl. te S t Disper- W. Plat S Disper- W. 

41595 36.0 1.14 

Mean equivalent width Mean equivalent width 


spectral type halfway between that Lacertae and Scorpii will 
about 0.5 equiv. The observed equivalent width Orionis, 09.5 from 
Table III, 0.28 equiv. Thus seems likely that the observed equivalent 
width real stellar atmosphere corresponding our O9.5 model 
atmosphere about 0.4 equiv. fairly certain that the equivalent width 
will not much larger than 0.5 equiv. much smaller than about 0.3 
equiv. When enter Fig. with these values, see that the ratio 
must greater than and probably lies between about and 25. That is, 
the present result indicates helium only 1/3 1/2 abundant pre- 
vious estimates (15), Aller (1), Bowen and Wyse (4), and Aller and 
Menzel (2) have indicated, and the reports the Transactions the 
International Astronomical Union (8) have suggested. 


| | 
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and the computed profile the model atmosphere, 13.3, are 
shown Fig. The intensity tracing for Lacertae from two-prism 
plate while that for Orionis from one-prism plate, mag- 
nified about make the dispersion the diagram the same for both 
tracings. mean profile from two more tracings drawn each case, 
although the grain from one plate. The shallow, extended Stark wings the 
computed profile can clearly traced the profile Lacertae 
where the line strong, and their existence suggested the profile for 
Orionis where the line weak and where the grain size and low dispersion 
make difficult obtain meaningful line profile. obvious that the 
computed line too strong for spectral type Thus, either have 


and the computed line profile when 13.8. 


assumed too much helium the model atmosphere, the conclusion reached 
above, wrong: spectral type has been assigned the model atmosphere, 
for some reason inherent the assumptions and calculations, the model 
atmosphere does not closely represent real atmosphere. 


order check the spectral classification the model atmosphere, the 
equivalent widths the CIII, and SilV lines used classify the model 
atmosphere (11) have been measured. Values are found essential agreement 
with those Aller (1), which were used the previous spectral classification 
(11) the model atmosphere. There little reason believe that the model 
atmosphere early Lacertae spectral type. Rather, halfway 
between Lacertae and Scorpii, that is, about Thus, since the 
spectral classification the model atmosphere appears correct, likely 
that have assumed too large abundance helium, and that the true 


There remains the question, only mentioned above, whether the structure 
model with small amount helium only will closely the same that for 
model with 5.67, and whether, indeed, the model represents real 
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stellar atmosphere all. order estimate the change structure caused 
different admixture helium, have integrated, but not tested, for 
radiative equilibrium, pure hydrogen model with the same temperature and 
surface gravity the model atmosphere. abbreviated summary 
the results given Table Both models have boundary temperature 
25,200° and log 4.200. assumed that the gas pressure given 
2P., and that the temperature law 


seen that the structures the two models not differ greatly, the 
electron pressure being reduced a_little the pure hydrogen model. The 
lowering the electron pressure with the lessening the helium content will 
tend advance the second ionization helium beyond the values used; how- 
ever, the electric fields arising from protons that cause the Stark effect the 
Hell lines increase slightly model with lessened helium content because 
the increased abundance hydrogen, hence there should little change 


TABLE 
THE DEPENDENCE LOG UPON HELIUM CONTENT THE O9.5 MODEL ATMOSPHERE* 


1.820 1.00 1.00 

0.05| 2.364 2.364 1.00 3.276 3.318 1.10 
2.513 2.513 1.00 2.20) 3.474 3.517 1.10 
2.911 1.05 3.606 3.657 1.12 


the the line-absorption coefficient lines Hell. conclude 
that the approximation using the structure model atmosphere computed 
with the ratio 5.67 compute the strength \4541 when 
the abundance helium rather less than this, valid and will not introduce 
serious errors into the theoretical curve growth. have confidence that the 
model atmosphere with 25,200° and log 4.200 gives good repre- 
sentation real O9.5 stellar atmosphere, since the strengths the CII, 
and SilV lines computed for spectral classification (11) are all 
the expected order magnitude. Since these lines have variety excitation 
potentials, and since the ions concerned have wide variety ionization 
potentials, seems unlikely that all the lines would have computed equivalent 
widths close the observed equivalent widths the structure the model 
atmosphere did not give fair representation real stellar atmosphere. For 
this model atmosphere, the effect the radiation pressure the mechanical 


| 
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equilibrium small; moreover the occurrence sharp-lined late O-type stars 
such Lacertae suggests that main-sequence late O-type stars have 
atmospheres mechanical equilibrium has been assumed. Thus, since the 
computations the strength the Hell line \4541 can brought into 
accordance with observation postulating reduced abundance helium, 
and since there are indications that the model not fairly representative 
real stellar atmosphere, conclude that the model atmosphere compu- 
tations presented here are evidence favor the ratio being between 
and about 25. That is, helium present only amount 
number, the rest the stellar atmosphere being effectively pure hydrogen. 


The recognition that the lines the Pickering series have shallow, 
greatly extended Stark wings, Fig. has some observational consequences. 
When visual inspection spectrograms made, will seem that stars 
with poor Jines, the lines, and particular, will tend 
disappear earlier spectral types than they will late O-type stars with 
good lines, for the eye more influenced the core line, than wide- 
spread, shallow wings. Thus since visual estimate the relative strength 
the lines used criterion spectral classification, there may 
systematic error the classification O-type stars with poor lines the 
effect that they are classified somewhat later than they should be. Also the 
shallow widespread wings the Pickering lines explain why these lines are 
difficult measure for radial velocity. 


For the model atmosphere (12) have found that the equivalent width 
13.33. Both these values are significantly less than the corresponding values 
for the model atmosphere, result that contradiction the observed 
relative strengths the Hell lines O-type stellar atmospheres. Observa- 
tionally has been established (see for instance Ref. (9)) that the Hell lines 
not reach maximum intensity until least and that the line \4541 
increases strength quite considerably from type O9.5 type O5. Either 
the abundance helium must greater early O-type stars than late 
stars, our model atmosphere does not give good representation 
pure absorption atmosphere. 


According Petrie (9) the equivalent width should expected 
least 1.0 equiv. O5. curve growth constructed with the 
equivalent widths \4541 computed above, Fig. and this diagram 
entered the point 20, which was found above fair estimate 
the model atmosphere, value that much too small. The relative strength 
\4541 the model atmosphere places the model about 
according Petrie’s (9) classification ratios. Such spectral type for the 
model consistent with the temperatures estimated for the early O-type stars 
from observation. The question then is, why are the calculated ionized helium 
lines weak for the model atmosphere? 
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will recalled that this model atmosphere (12) was assumed 
strict mechanical equilibrium and have 5.67 number. esti- 
mate the effect the structure the model altering the helium abundance, 
have integrated pure hydrogen model with the same temperature law 
and surface gravity. comparison the resulting electron pressure the 
the structure this model quite sensitive the assumed abundance 
helium and that the electron pressure sensibly lowered the abundance 
helium decreased. However, doubtful whether the changes due 
differing helium abundance are large enough account for the discrepancy 
found. seems more likely that the absorption-line early stars have atmos- 


TABLE 


= — | = = 


1.831 1.954 1.33 3.505 1.20 
0.10) 2.854 2.939 1.22 3.826 1.16 
2.966 3.049 1.21 3.907 3.968 1.15 
3.272 1.20 3.980 4.038 1.14 


found the appropriate function 


pheres state turbulence. Then the equivalent widths the absorption 
lines may larger (7, 18) than those predicted from stationary atmosphere 
strict mechanical equilibrium. When the computations the model 
atmosphere (12) are examined, seen that the radiation pressure exceeds 
the gas pressure through most the atmosphere. seems likely that under 
these circumstances strict mechanical equilibrium, such was forced the 
model, may not exist, and that the model atmosphere does not give good 
representation early O-type stellar atmosphere. Thus the computed 
strength the line the model atmosphere cannot safely 
used obtain estimate the relative abundance helium hydrogen 
early type stars. Our conclusions should based the results for the O9.5 
model atmosphere only. From the study this model atmosphere have 
concluded that the relative abundance number hydrogen helium lies 
the range 25. 
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